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1.  RESEARCH  OBJECTIVES  AND  SUMMARY 
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The  results  described  in  this  report  have  been  obtained  from  work 
whose  overall  objective  is  the  analysis  of  reacting  multiphase  flows  in 
advanced  air-breathing  propulsion  systems.  This  work  involves  the  formulation 
and  evaluation  of  new  theoretical  descriptions  of  multiphase  turbulent  mixing 
and  chemical  kinetic  phenomena  in  subsonic  and  supersonic  air  flows,  and  the 
generation  of  new  experimental  data  on  these  phenomena.-.  The  specific  object¬ 
ives  of  this  research  may  be  described  as  follows:  -i ,  ,<  - 

(a)  Spray  Modeling^  This  task  area  involves  the  development 
of  analytical  models  of  the  phenomena  occurring  in  the 
multiphase  turbulent  reacting  flows  involved  in  spray 
flame  evolution.  The  model  is  to  include  the  effects 
of  hetergeneous  and  homogeneous  finite  rate  chemical 
reaction,  turbulent  transport,  and  velocity  and  tempera¬ 
ture  nonequilibrium  between  the  droplets  and  gas  phase. 

Both  liquid  fuels  and  slurries  are  to  be  considered.  The 
analytical  modeling  work  is  to  be  supplemented  with  experi¬ 
mental  studies  of  spray  injector  characteristics  and  spray 
flame  structure  for  both  ramjet  and  gas  turbine  applications. 

(bj~* Modeling  of  Boron  Combustion  Processes^,  In  this  task  area 
advanced  analytical  models  of  the  combustion  of  boron  par¬ 
ticles  are  to  be  developed.  A  survey  of  the  existing  boron 
kinetics  literature  is  to  be  performed  to  provide  the  data 
base  for  the  development  of  advanced  boron  kinetics  models. 

This  survey  leads  to  the  definition  of  an  advanced  model  for 
boron  particle  combustion,  to  be  incorporated  in  modular 
models  of  slurry  combustion  process  in  ramjet.  Required 

boron  kinetics  data  are  also  to  be  defined.  -  „  ,  -  ) 

, jj  —  >  '  J 

(c)  .Modeling  of  Ramjet  Combustion  Processes^  The  modular  mod¬ 
eling  approach  to  the  analytical  characterization  of  high 
speed  combustion  processes  is  to  be  further  refined  and 
developed,  and  applied  to  the  analysis  of  results  genera¬ 
ted  from  experiments  relevant  to  high  speed  combustion 


processes.  Model  refinements  including  the  advanced 
quasiglobal  chemical  kinetics  model  for  the  character¬ 
ization  of  fuel-rich,  combustion  of  complex  hydrocarbon 
fuels  and  the  spray  model  formulations  developed  under 
item  (a)  above  are  to  be  incorporated.  This  modeling 
work  will  also  emphasize  slurry  fuels,  particularly 
boron  slurries.  Of  specific  importance  in  this  area 
is  the  experimental  data  on  gas-phase,  liquid-phase, 
and  slurry-fuel  combustion  phenomena  being  obtained 
at  the  Naval  Weapons  Center  under  a  joint  AFOSR/NWC 
sponsored  program.  Also  of  importance  in  this  area 
|  is  the  continued  development  of  unified  elliptic  com¬ 
putational  techniques  for  modeling  ramjet  combustion 
,  processes.  x" 

(d)  Modeling  of  the  Ducted  Rocket  Combustion  Processr  This 
task  area  involves  the  formulation  and  development  of 
a  modular  model  for  the  ducted  rocket  combustor.  This 
model  is  to  be  based  on  a  delineation  of  flow  charac¬ 
teristic  regions  obtained  through  the  use  of  three- 
dimensional  elliptic  aerodynamic  caculations,  and  is  to 
be  used  to  develop  a  broad  base  of  new  and  fundamental 
information  on  three-dimensional  multi-stream  interac¬ 
tion  in  reacting  ducted  flows. 


Work  carried  out  to  meet  these  objectives  during  the  contract  period  re 
suited  in  a  number  of  significant  accomplishments.  In  summary,  these  accom¬ 
plishments  include: 


(a)  Development  of  a  Comprehensive  Spray  Flame  Model:  A 
spray  flame  model  was  completed  which  includes  sub¬ 
models  for  two-phase  transport  and  diffusion,  hetero¬ 
geneous  and  homogeneous  combustion  kinetics,  turbulent 
mixing,  gas-phase  radiation,  and  buoyancy  effects. 
Based  on  the  assumption  of  near  dynamic  equilibrium, 
this  approach  allows  a  detailed  computation  of  the 
evolution  of  the  droplet  size  distribution  throughout 
the  flame,  and  incorporates  the  advanced  quasiglobal 
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kinetics  model  to  accurately  account  for  gas-phase 
hydrocarbon  combustion  processes  in  fuel-rich  regions 
of  the  flame.  Utilization  of  the  model  to  predict 
spray  flame  development  has  shown  the  necessity  of 
incorporating  accurate  representations  of  physical 
processes  which  control  the  flame  development,  in¬ 
cluding  those  mechanisms  which  occur  on  a  droplet 
scale  but  also  including  gas-phase  phenomena  invol¬ 
ving  finite-rate  kinetics,  thermal  radiation,  and, 
for  the  flowfield  considered,  buoyancy.  The  compu¬ 
ter  model  that  has  been  developed  from  this  work 
provides  a  generalized  framework  for  the  development 
and  testing,  in  conjunction  with  key  experimental 
data,  of  physically  perceptive  models  of  the  many 
mechanisms  involved  in  spray  combustion. 

(b)  Development  of  an  Accurate  Boron  Combustion  Model; 
Based  on  the  results  of  a  survey  of  boron  combustion 
modeling,  an  exiting  model  in  which  both  diffusion 
and  chemical  kinetics  mechanisms  are  implemented  was 
selected  for  further  development.  The  model  is 
based  on  a  clean-particle  assumption  -  it  is  assumed 
that  no  B2O.J  layer  is  present.  Since  particle  con¬ 
sumption  at  a  significant  rate  can  only  occur  for 
clean  particles,  this  approach  allows  examination 
of  the  effects  of  kinetics  limitations  on  boron 
consumption  rate.  The  model,  coupling  finite  rate 
oxidation  kinetics  and  diffusional  mechanisms,  shows 
excellent  agreement  with  experimental  data  for  boron 
particle  burn  time  as  a  function  of  particle  size, 
temperature  and  pressure.  Further,  the  model  indi¬ 
cates  that  kinetics  effects  enter  over  the  entire 
size  range  of  interest.  The  model  results  also  pre¬ 
dict  quenching  for  certain  combinations  of  particle 
size,  environment  pressure  and  temperature,  and 
environment  0^  concentration.  These  results  indicate 
a  critical  balance  between  heat  generation  and  heat 
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loss  <?n  the  particle  scale,  indicating  that  the 
coupling  of  fluid  dynamics,  and  particle  cloud  pro¬ 
cesses.  is  of  fundamental  importance  in  the  combus¬ 
tion  of  boron  in  practical  systems.. 

(c)  Modeling  of  Ramjet  Combustion  Processes;  Under  this 
program,  a  detailed  modular  model  for  sudden-expan¬ 
sion  ramjet  combustor  performance  evaluation  was 
completed.  This  model  combines  state-of-the-art 
representations  of  the  turbulent  mixing  and  chemical 
kinetics  rate  processes  with  submodels  for  the  fuel 
injection,  vaporization  and  mixing  phenomena  in  the 
combustor  inlet.  It  is  founded  on  the  use  of  a 
perfectly-stirred  reactor  model  to  represent  large- 
scale  recirculation  regions  in  the  combustor  and, 
in  particular,  the  flame  stabilization  characteris¬ 
tics  of  these  regions.  A  variety  of  computations 
have  been  carried  out  using  this  approach,  under 
this  and  other  Air  Force-  and  Navy-  sponsored  pro¬ 
grams,  to  analyze  and  interpret  experimental  data 
obtained  using  a  variety  of  sudden-expansion  ramjet 
configurations.  Results  of  this  work  have  demon¬ 
strated  the  role  of  the  sudden-expansion  recircula¬ 
tion  regions  in  flame  stabilization  at  lower  Mach 
numbers,  that  performance  fall -off  at  high  overall 
fuel-air  ratios  at  moderate  flight  Mach  numbers  can 
be  directly  related  to  the  effects  of  fuel-rich 
combustion  in  large  regions  of  the  combustor  flow- 
field,  and  that  a  desirable  fuel  distribution  at 
high  Mach  number  is  considerably  different  from  that 
which  is  desirable  at  lower  Mach  numbers.  All  of 
these  results  provide  a  critical  input  to  the  ana¬ 
lysis  and  interpretation  of  ramjet  performance  data. 

In  addition  to  the  modular  approach,  develop¬ 
ment  of  elliptic  Navier-Stokes  solution  techniques 
has  also  continued.  This  work  has  centered  on  the 
examination  of  turbulence  models,  and  has  shown  that 


the  algebraic  stress  modeling  approach  can  provide  a 
potentially  more  valid  description  of  turbulent  mixing 
in  combustor  flowfields.  The  increased  generality  of 
the  algebraic  stress  modeling  approach  is  of  particu¬ 
lar  importance  in  swirling  flows  where  large-scale  non¬ 
isotropy  of  the  apparent  viscosity  is  evident.  At  the 
same  time,  this  work  has  demonstrated  that  for  many 
combustors  of  interest,  for  example,  axisymmetric  sudden 
expansions  with  area  ratios  of  four  or  more,  inertia 
effects  dominate  in  much  of  the  flow,  thus  reducing  the 
effect  of  choice  of  turbulence  model  on  the  overall  flow- 
field  prediction. 

(d)  Ducted  Rocket  Model  Development:  Work  in  this  area  fo¬ 
cused  on  the  examination  of  different  techniques  for  the 
prediction  of  ducted  rocket  flowfields  and  performance. 

Aerodynamic  flowfield  predictions  were  carried  out  for  a 
representative  combustor  flowfield,  but  without  simula¬ 
tion  of  the  gas  generator  exhaust.  Limited  modular 
model  evaluation  of  ducted  rocket  performance  charac¬ 
teristics  was  also  accomplished.  Examination  of  avail¬ 
able  modeling  approaches  applicable  to  the  ducted  rocket 
showed  that  none  examined  accounted  for  the  mixing  and 
shock  interaction  effects  of  a  highly  underexpanded  gas 
gnerator  fuel  jet.  The  existence  of  this  key  feature  of 
the  flowfield,  which  may  well  dominate  the  overall  flow 
structure,  should  be  included  in  any  ducted  rocket  model- 
1 ing  approach. 

Although  as  noted  in  the  proceeding  paragraphs,  there  have  been  consider¬ 
able  advances  made  in  each  of  the  research  areas  considered  under  this  project, 
profitable  research  work  can  still  be  carried  out  in  each  of  these  areas.  In 
the  following  paragraphs  the  research  recommendations  that  result  from  the 
work  already  done  are  outlined: 

(a)  Spray  Modeling:  Two  avenues  of  approach  can  be  taken  in 
further  spray  modeling  work.  With  the  existing  model. 


experiments  can  be  designed  which  will  closely  approximate 
the  model  assumption,  so  that  direct  comparisons  of  the 
various  aspects  of  the  model  predictions  which  data  can  be 
made.  In  this  way,  each  of  the  submodels  incorporated  in 
the  overall  spray  model  can  be  examined,  to  provide  infor¬ 
mation  on,  for  example,  the  effects  and  importance  of  drop¬ 
let  interactions,  droplet/turbulence,  combustion/turbulence 
and  radiation/turbulence  interactions,  and  the  effects  of 
different  diffusion  rates  for  droplet  and  gas  phase.  This 
would  then  provide  a  wealth  of  further  data  on  droplet 
transport  and  spray  combustion  phenomena  in  a  regime  where 
the  droplet  and  gas  phase  motions  are  nearly  in  equilibrium. 
The  second  avenue  of  investigation  involves  the  relaxation 
of  the  near  equilibrium  hypothesis.  Within  the  context  of 
the  present  model  an  Eulerian  formulation  for  droplet  trans¬ 
port  can  be  constructed,  and  a  Lagrangian  formulation  can 
also  be  defined.  Note  that  by  carrying  out  such  an  analysis 
for  a  flowfield  which  in  general  satisfies  the  boundary 
layer  assumptions,  as  in  the  present  spray  flame  formulation, 
many  of  the  numerical  problems  inherent  in  elliptic  Navier- 
Stokes  solution  techniques  can  be  avoided,  an  important 
consideration  in  establishing  an  understanding  of  the 
phenomena  while  maintaining  relevance  to  airbreathing 
propulsion  systems. 

(b)  Boron  Combustion  Modeling:  Work  accomplished  to  date  shows 
the  importance  of  incorporating  both  diffusional  and  kinetics 
mechanisms  in  developing  an  accurate  model  for  particle-scale 
boron  combustion  processes.  There  are  several  areas  of 
research  and  model  improvement  that  this  work  suggest.  Ex¬ 
tension  of  the  model  to  include  the  effects  of  the  exist¬ 
ence  of  an  oxide  layer  is  clearly  a  critical  need  with 
respect  to  the  practical  use  of  the  existing  approach.  At 
the  same  time,  there  are  several  other  areas  which  require 
continued  research.  For  example,  the  potential  importance 
of  suboxide  formation  and  disposal  is  at  the  present  un¬ 
clear  and  requires  further  investigation,  while  more 
mechanistic  models  of  B-0,  B-O-H,  and  B-O-H-C  kinetics 


are  also  needed.  Further,  finite  rate  condensation  of 

B2O3  is.  important  in  the  overall  heat  release  process, 

and  itself  represents  an  important  area  of  continued 

* 

research  related  to  problems  of  boron  combustion. 

(c)  Modeling  of  Ramjet  Combustion  Processes:  Work  carried 
out  under  this  and  other  programs,  has  shown  the  funda¬ 
mental  advantages,  that  can  be  gained  through  the  use 
of  a  detailed  and  mechanistic  modular  model  of  the 
ramjet  combustion  process.  The  modular  model  developed 
in  this  work  is  however,  limited  to  nonswirling  sudden 
expansion  flowfields.  Similar  approaches  can  be  deve¬ 
loped  for  swirl  flow  in  subsonic  combustors  and  for 
supersonic  combustion  applications.  At  the  same  time, 
fundamental  model  development  work  is  required  using 
unified  Navier-Stokes  approaches,  and  expanding  the 
capabilities  of  these  models  to  characterize  flowfields 
with  chemical  reactions.  This  in  turn  indicates  a  need 
for  the  continued  development  of  simple  yet  sufficiently 
general  hydrocarbon  chemical  kinetics  models  as  well  as 
broadening  the  delineation  of  the  regimes  of  model 
applicability. 

(d)  Ducted  Rocket  Model  Development:  There  is  a  dearth  of 
information  on  the  mixing  characteristics  of  turbulent 
recirculating  flows  with  embedded  supersonic,  under¬ 
expanded  jets,  yet  this  type  of  combined  flowfield  is 
central  to  the  ducted  rocket  configuration.  Clearly 

a  need  exists  to  explore  the  characteristics  of  these 
flows,  both  experimentally  and  through  the  development 
of  theoretical  models,  in  order  to  develop  the  under¬ 
standing  necessary  to  control  them  in  practical 
applications. 


2.  RESEARCH  STATUS 


2.1  SPRAY  FLAME  MODEL 

The  spray  flame  model  developed  under  this  program  applies  to 
the  near-dynamic  equilibrium,  thermal  nonequilibrium  limit  of  a  two-phase 
flow.  Thus,  it  is  assumed  that  the  mean  velocity  of  both  the  droplets  and 
the  gas  phase  is  nearly  the  same,  but  diffusion  of  particles  relative  to 
the  gas  phase  is  allowed.  Gradient  diffusion  is  assumed,  with  an  apparent 
eddy  viscosity  that  can  be  different  from  that  which  applies  to  the  gas 
phase.  The  conditions  required  for  the  near-dynamic  equilibrium  assumption 
to  be  valid  are  clearly  not  met  near  the  exit  of  a  spray  nozzle  injecting 
fuel  into  a  quiescent  ambient.  Nevertheless,  the  extent  of  the  region  in 
which  this  assumption  is  not  valid  is  not  well  known,  and  it  can  be  expected 
to  become  more  valid  as  the  spray  f 1  arte  flowfield  develops.  Thus  it  is  of 
interest  to  investigate  the  capabilities  of  the  approach,  noting  that  when 
liquid  fuels  are  injected  using  air  streams  to  accomplish  primary  atom¬ 
ization,  the  near  dynamic  equilibrium  assumption  may  indeed  be  reasonably 
valid  over  most  of  the  flow. 

2.1.1  Governing  Equations 

Governing  equations  for  the  two-phase  spray  flame  flowfield  are 
obtained  following  the  development  outlined  in  Ref.  1.  A  key  element  of 
this  analysis  involves  the  definition  of  the  turbulent  fluxes  and  their 
relationship  to  the  gradients  of  the  appropriate  mean  quantities.  The 
"near"  equilibrum  concept  simply  recognizes  that  an  exact  equality  in 
velocities  between  the  phases  does  not  exist,  but  that  it  is  meaningful  to 
characterize  the  global  system  with  a  single  barycentric  velocity  in 
direct  analogy  with  molecular  diffusion.  It  is  assumed  that  the  droplets 
form  a  dilute  suspension  in  the  gas  flow,  and  can  be  treated  using  con¬ 
tinuum  concepts.  Thus  the  continuous  distribution  of  droplet  sizes  within 
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the  spray  flame  is  discretized  into  a  number  of  classes,  each  of  which 
contains  all  droplets  in  a  giver,  size  range  about  the  average  range 
for  the  class.  These  classes  are  theh  treated  as  droplet-phase  "species", 
and  transport  equations  are  written  for  each  class.  For  example,  if  the 
droplet  distribution  function  in  terms  of  droplet  radius  is  discretized 
into  a  number  of  radius  intervals,  the  independent  variable  becomes  the 
droplet  number  density  for  a  given  size,  defined  as  the  total  droplet  mass 
of  a  given  size  range  per  unit  volume.  The  total  droplet  density  is  ob¬ 
tained  by  summing  over  the  droplet  size  distribution  function,  and  the 
partial  and  total  particle  densities  are  determined  through  solution  of 
appropriate  transport  equations. 

The  diffusivities  of  the  droplets  will  generally  be  different 
from  that  which  pertains  to  the  gas  phase  components,  precluding  a  straight¬ 
forward  application  of  Fick's  Law. to  the  global  droplet-gas  system.  A 

multicomponent  diffusion  law  can  be  developed,  however,  by  treating  each 

> 

droplet  class  and  the  gas  phase  as  a  binary  subsystem  to  which  Fick's  Law 
may  be  appl ied. 

The  basic  governing  equations  are  then  as  follows;  where  the 
nomenclature  is  as  described  in  Table  3  following  Eq.  (17). 

Global  Continuity 
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Mixture  Energy  Equation  (cont'd). 
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Droplet  Energy  Equation  for  Droplet  Class 
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The  mixture  density  is  given  by  the  quation  of  state 

P  =  P/[Rot0^  ^ 
y  i  i 

where  is  the  molecular  weight  of  gas  phase  species  i,  Ro  is  the 
universal  gas  constant,  and  Tg  is  the  local  gas  temperature. 

To  close  the  equation  set  and  complete  the  model,  a  variety  of 

additional  submodels  are  required.  These  include  models  for  the  various 

production  rate  terms  w^-.w^j,  and  w^;  an  expression  for  the  thermal 

energy  transferred  to  the  droplets,  q  .  and  the  radiative  energy  transfer 
•  R 

from  the  system,  q  ;  a  droplet  tracking  scheme;  and  a  turbulence  model  to 
define  the  effective  gas-phase  eddy  di ffusi vi ty, evg. 


(7) 


Turbul ence  Model 


Turbulent  transport,  including  the  effects  of  buoyancy,  is 
modeled  using  a  two-equation  formulation  along  with  a  temperature  fluctua¬ 
tion  equation.  This  results  in  a  "three-equation"  formulation,  developed 
under  a  modified  version  of  the  Boussinesq  hypothesis  in  which  density 
fluctuations  are  neglected  everywhere  except  when  they  appear  in  a  buoy¬ 
ancy  term.  Writing  density  fluctuations  in  terms  of  temperature  fluctua¬ 
tions  leads  to  a  formulation  in  which  the  buoyancy  term  in  the  turbulent 
kinetic  energy  balance  is  related  to  the  root  mean  square  of  the  tempera¬ 
ture  fluctuation,  (Ref.  2)  Thus  the  t'.rbulence  model  is  given  Dy 
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Turbulent  Kinetic  Energy  Dissipation  Rate 
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and  the  temperature  fluctuation,  £ (T 1 )  ]  '  =  q  is  obtained  through  the 


solution  of  a 


Temperature  Fluctuation  Transport  Equation 


pu  |a  ,  pv  la  =  2 


_  1  3  /_  3q»  p  /3T» 

7F  '  7  '3?  (r  ^  tr‘  CglyT(3r' 


V  p7  * 


01) 


This  formulation  introduces  seven  adjustable  constants  and  three 
turbulent  exchange  coefficients.  For  most  of  the  work  described  in  this 
paper  a  single  set  of  constants  and  coefficients  was  used  in  the  turbulence 
model,  as  noted  in  Table  1.  The  basic  turbulence  model  coefficients:  Cp, 
C-j ,  C^,  and  ae  were  obtained  from  an  extensive  study  of  turbulent 
flowfield  modeling  reported  in  Ref.  3  while  the  coefficients  used  in  the 
buoyancy  modeling  C^,  C^,  C^,  C^,  and  aq  were  evaluated  through  compari¬ 
son  of  nonreacting  hot  jet  predictions  with  experimental  data. 


Table  1  Turbulence  Model  Constants 


Droplet  Tracking  Methodology 


The  spray  flame  model  requires  a  droplet  tracking  and  category 
readjustment  methodology,  since  the  rate  processes  result  in  droplet  size 
changes  and  these,  and  the  diffusion  of  droplets  of  different  classes, 
alter  the  droplet  distribution  at  each  local  point  in  the  flow.  At  each 
step  in  the  computation  a  new  average  size  for  each  class  must  be  defined, 
based  on  the  alterations  that  result  from  diffusion  and  rate  processes. 
This  procedure  follows  the  methodology  outlined  in  Ref.  4;  in  it  the 
appropriate  drop-size  average  is  defined  through  conservation  of  mass  and 
conservation  of  the  parameter  defining  the  dominant  rate  process  affecting 
droplet  size.  Thus,  if  the  rate  process  is  dependent  on  radius  to  the 
power  q,  the  appropriate  expression  for  the  new  average  radius  is  given  by 
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where  the  averaging  has  been  done  assuming  a  three-point  explicit  finite 
difference  scheme.  The  term  is  the  mass  fraction  of  droplet  class  j 

at  grid  point  n,m,  the  three  diffusional  components  that  contribute  to 
the  droplet  class  j  at  n+1 ,  m  are  ASpj-| »  A0p j 2  and  ^pj3;  ^pj  is  the  con~ 
tribution  from  finite  rate  processes.  Knowing  the  mass  fraction  and  size, 
the  number  density  of  droplets  of  class  j  at  point  n+1 ,m  is  then  determined 
from 
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Note  that  the  commonly  used  Sauter  mean  diameter  (SMD),  given  by  the 

3  2  2 

expression  d.^EN  -ds  /EN-d-  implies  a  d.  or  surface  area  dependent 

“TT  J  J  J  J  ^ 

rate  process,  i.e.,  q=2  in  Eq.  (12).  However,  the  ‘‘d  "  law  for  vapor¬ 
izing  droplets  implies  a  rate  process  that  is  proportional  to  the  droplet 
diameter,  so  that,  in  Eq.  (12)  q=l .  Somewhat  different  average  radii 
will  result  depending  on  the  governing  rate  process,  and  examination  of 
calculational  results  shows  that  the  use  of  q=2  (i.e.,  the  SMD  definition) 
in  Eq.  (12)  yields  a  slightly  different  overall  spray  development  than 
use  of  q=l  for  a  vaporizing  spray.  This  difference  is  small  for  the 
conditions  considered  in  this  report,  but  in  situations  in  which  several 
different  droplet-scale  rate  processes  are  occurring,  the  choice  of 
average  diameter  definition  may  become  more  critical.  Indeed  under  some 
circumstances  it  may  be  necessary  to  define  more  than  one  average  dia¬ 
meter  for  each  class  of  droplets. 


Thermal  Energy  Transfer  Terms 

Two  thermal  energy  transfer  terms  appear  in  the  mixture  energy 

•  R 

and  droplet-phase  energy  equations.  These  are  q  ,  appearing  in  the  mixture 
energy  equation,  which  represents  thermal  energy  transferred  out  of  the 
flame  by  radiation,  and  q  .  in  the  droplet  energy  equation,  which  repre- 
sents  a  transfer  of  thermal  energy  to  the  droplets  from  the  gas  phase. 
Radiation  from  the  flame  to  the  surroundings  is  assumed  to  be  dominated 
by  gas  phase  CO^  and  H^O  emission,  and  is  defined  by  an  empirical  expres¬ 
sion,  similar  to  that  used  in  Ref.  5: 


6.68x10  1?T4[aH  Qexp(-] .49xl0~3T)  + 
aCQ  exp{ -1 .81x10_3T)] 
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yielding  in  BTU/ft^-sec  with  T  in  °K. 


yielding  q^  in  BTU/ft^-sec  with  T  in  °K. 


Energy  transfer  between  gas  phase  and  particle  is  accounted  for 
through  relationships  obtained  (Ref.  6)  using  the  classical  spherical  droplet, 
flame  sheet  approach  and  solving  the  resulting  governing  equations  for 
the  droplet  heat  and  mass  transfer  processes.  For  a  vaporizing  droplet, 
and  accounting  for  droplet-gas  heat  transfer  resulting  from  both  conduc¬ 
tion  and  radiation,  this  yields  an  expression  for  the  energy  transfer 
to  the  droplet  of  the  form 
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Production  Rate  Terms 

Three  species  production  rate  terms  appear  in  the  formulation. 
These  involve  production  of  gas  phase  species  1  through  gas  phase  chemical 
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reaction,  w  -,  production  of  droplet  class  j  from  other  droplet  classes 
91  .  F  • 

through  the  effects  of  vaporization  and  combustion,  wp j ,  and  production 

of  gas  phase  species  i  through  vaporization  or  combustion  of  droplet  class 

j’  w 

The  mass  change  of  droplet  class  j  per  unit  volume  of  mixture  per  unit 
time  is  evaluated  using  the  classifical  spherical  droplet-flame  sheet 
model.  Thus  the  rate  of  change  of  droplet  mass  is  given  by 

*PJ  =  *n  (1  +  B)  06) 

P  L  j 

where  the  transfer  number,  B,  is  given  by 


B  =  r  [CJT  -LI  +  ffi&SJ 


The  species  produced  from  the  droplet  class  j  will  be  either  fuel,  in  the 
case  of  vapori  za  tion ,  or  CO-,  and  H^O  if  the  droplet  is  burning.  In  the 
model  described  in  this  report,  the  droplet  is  assumed  to  burn  if  suffi- 


cient  oxygen  is  present  in  its  immediate  environment;  otherwise  a  pure 
vaporization  process  takes  place.  Thus  if  otg  >0.001,  the  present 
model  produces  CO2,  1^0  (and  excess  fuel)  fro&  a  complete  combustion  pro¬ 
cess  in  the  flame  sheet,  and  otherwise  only  gas  phase  fuel  is  produced. 

This  lower  limit  is  of  course  arbitrary. 

Gas  phase  chemical  kinetics  are  evaluated  using  the  advanced 
quasiglobal  formulation  described  in  Ref.  7.  This  approach  includes 
models  for  the  fuel  pyrolysis  process  that  occurs  in  high  temperature, 
oxygen  deficient  regions,  and  thus  eliminates  the  inaccuracy  evident 
in  the  original  quasiglobal  model  formulation  (Ref.  8)  when  applied  to  fuel- 
rich  combustion.  In  the  advanced  quasiglobal  kinetics  approach,  the 
fuel  to  be  modeled  is  broken  down  into  aromatic  and  aliphatic  components, 
modeled  using  rates  developed  for  toluene  and  iso-octane,  respectively. 

That  is,  the  model  fuel  component  which  can  be  treated  as  aromatic  is 
assigned  rates  which  have  been  developed  by  extensive  comparison  of 
toluene  predictions  with  experimental  data,  and  the  model  fuel  component 
that  can  be  treated  as  aliphatic  is  assigned  the  rates  developed  for 
iso-octane.  Each  of  these  fuels  break  down  into  secondary  fuels,  C2H2 
from  toluene  and  C2H 4  from  iso-octane,  following  the  pure  pyrolysis 
pathway,  as  well  as  undergoing  partial  oxidation  by  attack  from  both 
02  and  OH.  The  secondary  fuels  then  react  with  the  remaining  OH  and 
O2,  and  the  wet  CO  mechanism  follows. 

For  the  computations  to  be  described  in  this  report,  the  fuel  was 
pentane,  so  that  only  the  aliphatic  component  needed  to  be  included  in 
the  kinetics  formulation.  Since  NOx-formation  reactions  were  also 
neglected,  the  gas-phase  kinetics  model  used  for  these  computations 
involved  the  20  reactions  among  the  10  species  shown  in  Table  2. 


TABLE  2.  ADVANCED  QUASI6L0BAL  KINETICS  MODEL  FOR 
ALIPHATIC  FUELS 
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(b)  ELEMENTARY  STEPS 
k,  --  ATBexp(-E/RT) 
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TABLE  3.  NOMENCLATURE  FOR  EQ'S.  1-17 


B  transfer  number,  defined  by  Eq.  I? 

Cp  liquid  phase  specific  heat 

C|  1 

Cj  I  coefficients  in  turbulent  kinetic  energy  dissi- 
C1  (  pation  rate  equation,  defined  in  Table  !. 

C.  > 

!  coefficients  in  temperature  fluctuation  equa- 
Cq  j  tion,  defined  in  Table  1. 

Cp  eddy  viscosity  coefficient,  defined  in  Table  1 
e  droplet  ennssivity 
g  gravitational  constant 

H  Mixture  total  enthalpy  per  unit  mass  of  mixture 
--  h  ♦  I  fi  h  .  ♦  u:/2 
9  j  PJ  PJ 

h  Static  enthalpy  of  gas  phase  per  unit  mass  of 
y  mixture  =  I  o  h  ( T) 
i  1  1 

h  Static  enthalpy  of  droplet  class  j  per  unit 
mass  of  droplet  class  j 
J  ■  gas  phase  species  i  diffusional  flux 

Jpj  droplet  phase  species  j  diffusional  flux 

k  turbulent  kinetic  energy 

kf  forward  reaction  rate  coefficient,  Table  2 

L  latent  heat  of  vaporization  per  unit  mass 
tT  number  density  of  droplets  of  class  j,  Eq.  13 

Nu  Nusselt  number  for  gas-droplet  heat  transfer 
P  pressure 

Pr  turbulent  Prandtl  number 

q  root-mean-square  temperature  fluctuation 

qr  energy  transfer  from  mixture  by  gas-phase 
radiation,  £q.  14 

qpj  energy  transfer  to  droplet  class  j,  Eq.  IS 

0  heat  of  vaporization 

r  radial  coordinate 

r j  radius  of  droplet  of  class  j 

Ro  universal  gas  constant 

Sc  turbulent  Schmidt  number 

Tg  gas  phase  temperature 

T  gas  phase  temperature  fluctuation 
T^  droplet  temperature 

u  axial  velocity 
v  radial  velocity 

*oi  production  of  gas  phase  species  i  through  chem- 
f  leal  reaction 

»  .  production  of  droplet  phase  class  j  through 
^  effects  of  vaporization,  combustion,  ar,d  agglo¬ 
meration,  per  unit  volume  of  mixture  per  unit 
p  time 

i  ■  ■  production  of  gas  phase  species  i  ti. rough 

J  vaporization  or  combustion  of  droplet  cl.-ss  j, 
per  unit  volume  of  mixture 
kh  molecular  weight  of  species  i 
«  axial  coordinate 

To  ctygen  concentrat ion,  fraction  of  cas  phase 
jg  mixture  mass  fraction  of  gas  phase  =  Eat 

i  mixture  mass  fraction  of  species  i 
3  mixture  mass  fraction  of  droplet  class  j 


ratio  of  droplet  class  j  and  gas  phase  thermal 
diffusivities  4  e„  r 

Tj  P9 


turbulent  kinetic  energy  dissipation  rate 
droplet  phase  thermal  diffusivity,  class  j 

gas  phase  eddy  viscosity 

gas  phase  thermal  conductivity 
turbulent  effL-tive  viscosity  =  pc 

mixture  density 
liquid  density 

ambient  density 
Stefan-Boltzmann  constant 

ratio  of  eddy  viscosity  of  droplet  class  j  to 
gas  phase  eddy  viscosity  =  r  /r 

VJ  v9 

effective  "Prandtl  number"  for  turbulent  kinet¬ 
ic  energy.  Table  1 

effective  "Prandtl  number"  for  temperature 
fluctuation.  Table  1 

effective  "Prandtl  number"  for  turbulent  energy 
dissipation.  Table  1 . 
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2.1.2 


Comparison  with  Experimental  Results 


Vaporizing  Spray 

In  Ref.  9,  Shearer,  et  al .  describe  detailed  measurements  of 
mean  velocity,  total  concentration  (both  liquid  and  vapor  phase),  and 
Reynolds  stress  which  were  made  in  a  two-phase  Freon-ll/air  jet  mixing 
with  still  air.  The  jet  was  produced  by  the  injection  of  Freon-11 
into  the  air  using  an  air-atomizing  nozzle;  the  initial  1 iquid-phase 
Freon  mass  fraction  was  0.87.  Shearer,  et  al .  do  not  report  droplet 
distributions;  however,  they  give  the  initial  SMD  (Sauter  mean  diameter) 
of  the  Freon-11  droplets  as  29  microns.  The  injection  temperature  of 
the  Freon-11  was  just  below  its  saturation  temperature,  which  in  turn 
was  just  below  the  ambient  air  temperature.  Thus  the  vaporization  pro¬ 
cess  for  this  jet  was  driven  primarily  by  concentration  gradients. 

Within  the  context  of  the  model,  this  vaporization  process  was  approxi¬ 
mated  by  assuming  a  transfer  number,  8.  of  ten,  in  the  droplet  vapori¬ 
zation  expression,  Eq.  (16).  Vaporization  was  assumed  to  initiate  when 
the  droplet  temperature  reached  the  saturation  temperature.  An  assumed 
initial  droplet  distribution,  obtained  from  literature  data  for  air- 
atomizing  nozzles,  was  used  to  initiate  the  calculations;  this  distribu¬ 
tion,  shown  in  Fig.  1,  has  a  29  p  SMD.  All  other  initial  conditions 
(i.e.,  velocity,  temperature  and  droplet  mass  fraction,  turbulence  energy 
and  dissipation  rate)  were  taken  to  be  uniform  across  the  jet.  All  drop¬ 
let  classes  were  taken  to  diffuse  at  the  same  rate  as  the  gas  phase,  i.e., 
6-=l  for  all  j.  The  phenomenon  of  phase  separation,  or  relative  diffusion, 

J 

is  included  in  the  formulation  but  in  the  absence  of  detailed  data  8j 
must  be  established  parametrical ly,  so  that  8-=l  provides  a  reasonable 

J 

baseline  case. 

Axial  variations  of  centerline  velocity  and  species  mass  fractions 
are  cor, pared  with  the  data  of  Shearer,  et  al .  in  Fig.  2.  As  noted  on 
the  figure,  these  computations  have  been  carried  out  for  a  turbulent 
Prandtl  number  of  0.70.  These  results  indicate  that  the  predicted  velo¬ 
city  profile  compares  reasonably  well  with  the  experimental  results. 


but  that  the  centerline  Freon  concentration  (both  liquid  and  vapor  phase) 
is  somewhat  underpredicted.  More  accurate  predictions  of  the  centerline 
Freon  concentration  can  be  achieved  by  using  a  unity  turbulent  Prandtl 
number  assumption,  but  results  obtained  with  this  assumption  strongly 
underpredict  the  width  of  the  concentration  profile.  As  can  be  seen 
from  Fig.  3,  the  present  computations  (at  a  turbulent  Prandtl  number  of 
0.70)  provide  very  accurate  predictions  of  the  shape  and  width  of  the  con¬ 
centration  profiles. 

Radial  profiles  of  mean  axial  velocity  are  also  predicted  very 
well  by  the  vaporizing  spray  jet  model,  as  shown  by  the  results  pre¬ 
sented  in  Fig.  4.  Note  that,  despite  its  two-phase  nature,  the  vaporizing 
spray  jet  exhibits  the  profile  similarity  behavior  expected  of  an  axisym- 
metric  single-phase  jet  into  still  air. 

The  analytical  model  described  here  provides,  in  addition  to  the  mean 
profile  results  already  discussed,  a  variety  of  information  with  respect  to 
the  development  of  the  spray.  For  example,  the  evolution  of  the  drop  size 
distribution  in  the  spray  is  an  output  of  the  computation,  and  results  for 
this  parameter  are  shown  in  Fig.  5.  The  drop  size  distribution  evolution 
is  shown  here  for  the  centerline  and  an  r/x  value  of  0.10,  which  corresponds 
roughly  to  the  half-radius  of  the  spray  jet.  Of  interest  is  the  fact  that 
as  the  number  of  droplets  decreases,  the  distribution  remains  roughly  the 
same,  so  that  the  SMD,  for  example,  shows  little  change.  The  nearly  verti¬ 
cal  lines  on  Fig.  5  represent  the  average  diameter  for  each  class  at  each 
value  of  X/D;  the  fact  that  they  are  nearly  vertical  indicates  little 
change  in  mean  diameter  is  occurring  despite  a  nearly  two-order-of-magni tude 
decrease  in  droplet  number  density.  Thus  the  mean  diameter  of  a  spray  dis¬ 
tribution  does  not,  from  these  results,  appear  to  be  a  good  indicator  of  the 
spray  vaporization  rate  by  itself. 

Burning  Spray 


Mao,  et  al.  (Ref.  10)  report  results  for  a  spray  diffusion  flame.  The 
experimental  apparatus  and  measurement  techniques  used  are  quite  similar  to 
those  used  in  the  vaporizing  spray  study  of  Shearer,  et  al  .  already  dis- 
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igure  3.  Radial  Variation  of  Freon  Concentration 
Vaporizing  Spray  Jet.  Data  from  Ref.  9 
Model  Predictions  at  X/D  =  235 


drop  diameier,  microns 


igure  5.  Computed  Evolution  of  Droplet  Size 
Distribution,  Vaporizing  Spray  Jet 


cussed.  For  the  spray  flame  experiments,  Mao  et  al .  utilized  n-pentane  as 
a  fuel,  with  an  initial  droplet  distribution  having  a  SMD  of  35  pm.  Measure¬ 
ments  include  centerline  velocity,  fuel  species,  and  temperature  distributions, 
as  well  as  radial  profiles  of  velocity  and  temperature.  Analysis  of  these 
data  using  the  1 ocally-homogeneous-flow  (LHF)  approximation  is  also  reported 
in  Ref.  10. 

The  calculations  :arried  out  using  the  spray  flame  model  described 
in  this  report  were  initiated  using  a  droplet  distribution  similar  to  that 
used  for  the  nonburning  spray,  but  adjusted  to  yield  an  SMD  of  35  pm.  Ini¬ 
tiation  of  the  flame  in  the  calculations  is  provided  by  the  droplet  flame 
sheet  assumption:  since  an  air-atomizing  nozzle  was  used  sufficient  oxygen 
is  available  at  the  injector  exit  to  initiate  chemical  reaction.  This  is  an 
artifice,  but  it  does  provide  an  appropriate  flame  initiation  zone  to  begin 
the  computations.  This  initial  reaction  quickly  exhausts  the  available  oxy¬ 
gen,  so  that  the  bulk  of  the  droplets  vaporize  in  a  hot  fuel  rich  region  in 
the  interior  of  the  flame.  The  fact  that  in  the  computations  the  bulk  of  the 
droplets  do  not  enter  directly  into  the  combustion  process  is  fully  supported 
by  experimental  observations  made  by  Mao  et  al .  Gas  phase  combustion 
is  described  in  the  model  using  the  advanced  quasiglobal  approach  described 
in  the  preceding  section,  and  gas  phase  thermal  radiation  is  also  included. 

This  latter  phenomenon  was  found  to  strongly  affect  the  development  of  the 
flame,  as  described  below. 

A  comparison  of  predicted  centerline  temperature  with  the  temperature 
experimentally  measured  in  the  pentane  spray  flame  is  shown  in  Fig.  6. 

The  overall  profile  shape  is  very  well  predicted,  particularly  for  X/D>200, 
but  the  predicted  temperature  is  everywhere  higher  than  that  measured. 

Here  radiation  may  play  a  role,  both  in  the  model  and  in  the  experiment. 

For  X/D<200,  the  predicted  temperature  is  substantially  in  excess  of  that 
measured.  In  these  calculations,  the  droplets  vaporized  by  X/D= 1 00 , 
whereas  experimental  observations  of  droplets  to  150  diameters  were  reported. 
This  suggests  that  the  computed  droplet  vaporization  rate  is  too  high,  and 
that  the  flame  sheet  model  for  droplet  consumption  on  the  edge  of  the  flame 
is  an  oversimplification.  Nevertheless,  noting  that  there  is  no  normaliza¬ 
tion  in  this  comparison  (other  than  division  of  the  axial  distance  by  the 
physical  jet  diameter),  these  results  are  very  encouraging. 


Both  thermal  radiation  and  buoyancy  effects  had  to  be  included  in 
the  computational  model  in  order  to  achieve  a  reasonable  representation  of  the 
spray  flame  development.  A  number  of  parametric  computations  were  carried 
out  during  the  development  of  the  approach,  the  results  of  which  illustrate 
the  effects  of  the  neglect  of  these  important  phenomena.  For  example,  the 
results  shown  in  Fig.  7  indicate  centerline  temperature  distributions  ob¬ 
tained  without  including  the  effects  of  buoyancy,  but  with  and  without  in¬ 
corporation  of  the  effects  of  gas-phase  thermal  radiation.  The  measured 
centerline  temperature  peaks  at  a  value  of  some  1800°K.  Without  incorpora¬ 
ting  the  effects  of  gas-phase  thermal  radiation,  the  predicted  centerline 
temperature  continues  to  increase,  eventually  reaching  the  adiadatic 
flame  temperature  for  stoichiometric  pent^ne-air  combustion.  Incorporating 
radiation  losses,  the  model  prediction  of  centerline  temperature  is  much 
reduced  with,  in  this  case,  a  predicted  temperature  peak  of  only  1550°K. 

Buoyancy  influences  both  the  mean  flow  and  turbulence  structure  in 
the  spray  flame,  and  the  relative  effects  of  buoyancy  increase  as  the  magni¬ 
tude  of  the  buoyancy-induced  body  force  term  in  the  momentum  equation  increases 
relative  to  the  momentum  flow  per  unit  area  term.  Since  in  a  jet-i nto-stil 1  - 
air  the  local  momentum  flow  per  unit  area  is  continuously  decreasing  with 
axial  distance,  at  the  same  time  as  the  temperature  and  therefore  body 
force  term  increases,  buoyancy  effects  can  be  expected  to  always  exert  some 
influence  on  the  development  of  a  spray  flame  in  a  quiescent  ambient.  The 
magnitude  of  this  influence  is  illustrated  by  the  results  shown  as  Fig.  8, 
which  presents  computed  centerline  temperatures  obtained  including  the 
effects  of  gas-phase  thermal  radiation,  but  with  and  without  incorporation 
of  the  effects  of  buoyancy  as  the  mean  flow.  Note  that  these  results 
do  not  include  the  influence  of  buoyancy  effects  on  the  turbulent  mixing 
rate.  The  results  shown  in  Fig.  8  show  that  there  are  very  large  effects 
of  buoyancy  in  this  spray  flame  structure,  since  with  mean-flow  buoyancy 
effects  included  the  calculation  represents  the  shape  of  the  centerline 
temperature  profile  quite  well,  particularly  with  respect  to  the  rapid  decay 
of  centerline  temperature  downstream  of  the  location  of  the  temperature. 
However,  the  effects  of  incorporation  of  mean  flow  buoyancy  without  incorpora¬ 
ting  the  corresponding  buoyancy-turbul ence  structure  effects  produce  a  quite 


different  result  when  computed  centerline  velocity  is  compared  with  experi¬ 
mental  data.  This  is  demonstrated  by  the  results  shown  in  Fig.  9  where  it 
can  be  seen  that  with  mean-flow  buoyancy,  an  increase  in  centerline  velocity 
is  observed  for  100<X/D<  300.  Clearly  the  data  do  not  support  this  result. 
However,  if  the  effects  of  buoyancy  on  the  turbulent  mixing  process  are 
included,  through  the  buoyancy  terms  appearing  in  Equations  (9)  and  (10), 
this  centerline  velocity  increase  phenomenon  is  to  a  large  extent  avoided. 
This  result  is  shown  by  the  centerline  velocity  variation  results  shown 
in  Fig.  10,  which  correspond  to  the  centerline  temperature  results  of 
Fig.  6.  These  results  show  that  the  effects  of  buoyancy  become  apparent  at 
around  X/D=200,  and  result  in  a  reduced  rate  of  centerline  velocity  decay 
compared  to  the  experimental  data.  The  underprediction  of  the  velocity 
decay  rate  is  consistent  with  the  overprediction  of  the  centerline 
temperature,  one  effect  of  which  would  .be  to  increase  the  magnitude  of 
the  buoyancy  effect.  Yet  velocity  profiles.  Fig.  11,  agree  extremely 
well  with  the  data  for  X/D>200;  these  profiles,  while  normalized  with 
the  centerline  velocity  (as  was  the  data)  are  not  normalized  with  respect 
to  the  physical  coordinate  in  a  way  that  would  mask  disagreement.  Indeed, 
the  normalization  used  indicates  that  the  computations  correctly  model 
the  spread  rate  of  the  spray  diffusion  flame,  at  least  for  X/D>170. 

The  accuracy  of  the  spread  rate  prediction  is  also  indirectly 
supported  by  comparison  of  the  calculated  and  measured  Reynolds  stress 
component.  This  result,  shown  in  Fig.  12,  shows  a  very  good  level  of 
agreement  both  in  magnitude  and  location  of  the  maximum  shear  stress.  It 
has  been  established  for  single  phase  jets  into  quiescent  ambients  that 
the  rate  of  spread  of  the  jet  is  controlled  by  the  peak  shear  stress  levels 
encountered  in  the  profile,  so  that  the  good  level  of  agreement  seen  in 
Fig.  12  supports  the  velocity  profile  agreement  displayed  in  Fig.  11. 

Finally,  computed  and  measured  temperature  profiles  are  compared 
in  Fig.  13.  Here  it  is  evident  that  in  the  center  of  the  flame  the  model 
predicts  higher  temperatures  than  are  experimentally  observed.  Again,  both 
droplet  vaporization  rate  and  gas  phase  radiation  submodels  are  implicated 
in  this  d i sagreement .  Other  factors  may  enter  as  well.  For  example, 
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Figure  12.  Reynolds  Stress  Profile,  Reacting 
n-Pentane  Jet 


although  the  flow  has  been  treated  as  fully  turbulent,  in  both  the  original 
experiment  documentation  and  in  this  work,  it  may  have  actually  involved 
a  laminar  or  transitional  flow  in  the  first  100  diameters  or  so.  If  this 
were  the  case,  a  slower  rate  of  mixing  than  predicted  herein  would  have 
existed,  resulting  in  lower  centerline  temperature  than  predicted  by  this 
fully  turbulent  model. 

2.1  .3  Pi scussion 

In  assessing  the  performance  of  the  spray  flame  model  just  des¬ 
cribed,  it  is  important  to  keep  in  mind'  the  important  effects  of  the  the 
strong  coupling  inherent  in  a  spray  diffusion  flame.  For  the  particular 
flow-field  considered  this  coupling  includes  not  only  the  mixing  and  kin¬ 
etics,  but  also  the  effects  of  phase  change,  radiation,  and  buoyancy. 

Taking  just  one  of  these  as  an  example,  the  effects  of  inclusion  of  thermal 
radiation  were  found  to  be  significant  in  these  calculations,  accounting 
for  centerline  temperature  differences  ranging  from  400°K  to  1000°K.  This 
level  of  temperature  difference  is  significant  with  respect  both  to  chemical 
kinetic  and  buoyancy  effects.  The  buoyancy  model  is  also  worth  noting, 
since  results  of  additional  computations  clearly  showed  that  neglecting 
the  influence  of  buoyancy  on  the  turbulent  structure  produced  erroneous 
results.  Thus,  while  attention  in  assessing  the  performance  of  this  model 
may  be  focused  on  the  use  of  the  near  dynamic  equilibrium  hypothesis,  it 
is  not  altogether  clear  that  this  hypothesis  exerts  a  dominant  effect  on 
the  overal 1  resul ts . 

Overall,  the  approach  described  in  this  paper  provides  a  reasonably 
valid  picture  of  the  development  of  a  spray  diffusion  flame.  Moreover,  it 
provides  a  computational  tool  which  can  usefully  be  applied  to  the  parametric 
examination  of  the  many  physical  and  chemical  interactions  that  underlie 
spray  flame  phenomena,  both  from  a  purely  theoretical  standpoint  and  in 
conjunction  with  detailed  spray  flame  measurements.  The  approach  also  pro¬ 
vides  the  basis  for  further  development,  for  example  through  relaxation  of 


the  limitation  to  near  dynamic  equilibrium  in  the  spray  transport  process. 
Clearly,  the  model  is  complex,  and  involves  a  great  variety  of  phenomena. 

But  spray  diffusion  flames  are  themselves  complex,  and  careful  theoretical 
and  experimental  study  is  required  to  determine  where  model  simplification 
may  justifiably  be  made. 

From  a  computational  standpoint  the  complexity  of  the  model  is 
actually  more  apparent  than  real.  Even  with  detailed  finite-rate  chemistry 
as  described  by  the  quasiglobal  kinetics  model,  computation  of  the  n-pentane 
spray  flame  considered  here,  carried  out  to  an  axial  distance  of  500  nozzle 
diameters,  requires  only  about  30  minutes  of  CPU  time  on  a  CDC  7600.  Com¬ 
putational  times  on  more  advanced  computers  can  be  expected  to  be  consider¬ 
ably  less.  What  the  model  provides  for  the  investment  in  computer  time  is 
substantial  detail  with  respect  to  the  general  structure  of  the  spray-flame 
flow  field.  For  example,  detailed  temperature  contours  in  the  spray  flame 
can  be  presented  to  illustrate  the  structure  of  the  flame,  as  in  Fig.  14. 

When  compared  with  dropl  et  di  stributions ,  as  shown  in  terms  of  total  mass 
fraction  in  liquid  form  on  Fig.  15,  conclusions  can  be  drawn  as  to  the 
mechanisms  which  the  modeling  indicates  will  dominate  in  the  flame  struc¬ 
ture.  Comparison  of  Figs.  14  and  15  shows  that  the  bulk  of  the  droplets 
are  found  in  regions  on  the  flame  in  which  local  temperature  is  greater 
than  800°K,  £>ut  below  1200°K.  Comparison  with  gas-phase  oxygen  profiles 

i 

shows  that  this  is  a  very  fuel-rich  region.  Thus,  except  for  an  initial 

short  region  near  the  nozzle  exit  where  temperatures  range  from  1600°K  to 

2000°K  (and  which  is  an  artifice  of  the  assumption  that  droplets  burn, 

rather  than  vaporize,  ifctn>0.001),  vaporization  is  the  primary  mechanism 

u2 

responsible  for  the  reduction  in  liquid-phase  mass  fraction.  Combustion 
occurs  only  in  a  narrow  region  surrounding  the  droplet  cloud,  and  for  X/D<80, 
is  a  low-temperature,  fuel -rich  process  involving  the  mechanisms  of  pure 
and  oxidative  pyrolysis.  Beyond  X/D-80,  the  liquid-phase  mass  fraction 
drops  below  0.05  and  there  is  a  sudden  and  rapid  increase  in  the  reaction 
rate  as  indicated  by  the  widening  of  the  predicted  flame  shape  seen  on 
Fig.  14.  This  behavior  is  qualitatively  in  agreement  with  that  observed 
experimentally,  although  these  predictions  indicate  its  occurrence  rather 
sooner  in  the  spray  flame  development  than  was  otherwise  observed. 
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It  should  also  be  noted  that  the  particular  experimental  con¬ 
figuration  selected  for  study  in  this  work  was  not  ideal  from  the  stand¬ 
point  of  modeling  complexity.  In  particular,  the  quiescent  ambient  intro¬ 
duces  larger  effects  of  buoyancy  than  would  be  the  case  in  a  coflowing 
situation.  Indeed,  it  is  possible  to  design  an  experiment  which  satisfies 
the  assumptions  of  the  analysis  much  more  closely  than  that  considered  here, 
for  example,  with  respect  to  the  near-dynamic  equilibrium  assumption,  while 
at  the  same  time  reducing  the  magnitude  of  the  buoyancy  effects.  By  doing 
this,  and  carefully  comparing  the  experimental  results  with  analytical  pre¬ 
dications  it  would  be  possible  to  considerably  enhance  the  characteristics 
of  the  multiphase  spray  flame  model  described  in  this  work.  By  inference, 
then,  questions  of  droplet-scale  modeling,  droplet-droplet  and  turbulence 
interaction  could  all  be  addressed  in  a  step-by-step  manner. 

Finally,  the  results  presented  here  are,  in  terms  of  level  of 
agreement  with  the  data,  substantially  similar  to  those  reported  in 
Ref.  10.  The  results  described  in  Ref.  10  were  obtained  using  a  model  which 
ignores  two-phase  flow  effects,  finite-rate  vaporization,  finite-rate  gas 
phase  kinetics,  gas-phase  thermal  radiation,  and  buoyancy  phenomena. 

Instead,  the  model  utilized  in  Ref.  10  incorporates  a  turbulence-chemistry 
interaction  model  based  on  an  assumed  probability  distribution  function 
approach.  Two-phase  flow,  finite-rate  vaporization,  finite-rate  gas  phase 
kinetics,  gas-phase  thermal  radiation,  and  buoyancy  are  all  physically  real 
phenomena,  and  the  parametric  investigation  reported  above  indicate  the  mag¬ 
nitude  of  their  effects.  Thus,  the  agreement  with  experiment  of  a  theoretical 
treatment  that  ignores  each  of  these  phenomena  in  succession  requires  further 
explanation. 

2.2  MODELING  OF  BORON  COMBUSTION  PROCESSES 

Boron  slurry  fuel  combustion  is  a  complex  phenomenon  involving 
two-phase  flow  and  finite-rate  reaction  processes  in  the  aas  phase  and  in 
the  particle  surface.  In  a  slurry  fuel,  the  boron  is  generally  in  the  form 
of  micron-sized  particles  suspended  in  a  carrier  hydrocarbon.  As  the  hydro¬ 
carbon  carrier  burns,  the  boron  is  heated  to  the  ignition  temperature,  and 


it  then  subsequently  burns  with  exposure  to  sufficient  oxygen.  One  of  the 
difficulties  encountered  in  boron  combustion  is  that  a  surface  oxide  layer 
can  form  on  the  particle.  This  oxide  has  a  boiling  point  below  that  of  the 
metal,  so  that  the  boron  combustion  process  involves  a  diffusion  of  oxidizing 
gases  through  the  liquid  oxide  layer  with  a  surface  reaction  taking  place 
at  the  particle/liquid  interface.  Since  this  is  a  slow  process,  one  strat¬ 
egy  for  obtaining  high  boron  combustion  efficiencies  is  to  limit  the  amount 
of  oxygen  reaching  the  surface  of  the  particle  during  the  heat-up  process; 
i.e.,  carrying  out  the  particle  heatup  in  a  regime  which  is  fuel -rich  rela¬ 
tive  to  the  hydrocarbon  carrier,  and  then  exposing  the  hot  particles  and 
remaining  fuel  to  additional  air.  This  form  of  staged  combustion  is  diffi¬ 
cult  to  achieve  in  practice,  and  the  development  of  boron  slurry  fueled 
combustors  could  be  materially  aided  by  the  availability  of  analytical  design 
tools.  The  modular  model  described  in  the  next  section  of  this  report  pro¬ 
vides  the  framework  for  the  type  of  design  tool  that  could  be  used  for  this 
purpose,  but  its  extension  to  the  analysis  of  boron  slurry  combustion  re¬ 
quires  the  development  of  submodels  of  the  boron  combustion  process.  These 
submodels  are  the  principal  focus  of  the  work  carried  out  under  this  portion 
of  the  AFOSR-sponsored  research  at  SAI. 

2.2.1  Boron  Combustion  Models 

Single-particle  boron  ignition  and  combustion  models  have  been 
reported  by  King  (Refs.  11,  12),  Meese  and  Skifstad  (Ref.  13),  Vovchuck  et 
al.  (Ref.  14),  Mohan  and  Williams  (Ref.  15)  and  Edelman,  et  al .  (Refs.  16, 

17).  The  component  equations  that  make  up  each  of  these  models  are  sum¬ 
marized  in  Table  3,  which  also  includes  a  review  of  the  assumptions  built 
into  each  of  the  models.  The  nomenclature  for  Table  3  is  given  in  Table  4. 

In  general,  these  formulations  are  based  on  the  assumption  of  a  spherical 
particle  (except  for  the  work  of  Mohan  and  Williams)  and  on  the  conserva¬ 
tion  equations  for  mass  and  energy.  The  momentum  equation  is  not  explicitly 
coupled  into  the  equation  set,  but  in  certain  of  the  models  the  relative 
velocity  effect  is  accounted  for  through  a  Nusselt  number  correction  on  heat 
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lOJAIJMS  FOR  h 


nw*-.;  hc-hs 


F  rq.  tnvKdrs:  an  Arrhenius  *  1  s  c  es  1 ty  - 1«  "pc  r  a  1  ur  c  law,  V  -  A  » 


Henry’s  law,  Fl/U  -  Kl 

WlUr  i  or  r  r  1*  1 1  on,  0|B  -  7.4  *  30"’ 


(4  ^  r  7 


l^JAl ) OK  MUM  NCtS 


o  Siren's  constant  -  5  6C9  •  10"*  U/oi*  K  •  0.17)4  •  )  0’  *  ElU/h-ft1  Kfc 


D-D0(T/I0)'  ”  X  -  >o  (l/7o)  — 


I,  ■  U  1  /  k  P  -  r.p  (H/R‘7.  -  H/R*T  ) 

O  t  DC 

K  D  can  be  expressed  In  an  Arrhenius  form  (sec  King) 


Clausius  -  Claptryon 


h  h 


h  -  0.347  -  3  0"1  (hu)  7^*/  (rp  *  X)  (1S74)  h  *  0 .194  •  10" 1  7^*/  r?  (J97Z) 

h  -  (1  .  0. ?9L  Prl/’  Re17’)  t^/r 

The  first  stage  of  a  two  stage  burn.  During  this  tlr-e  thp  oxide  layer  Is  removed. 

A  state  at  which  r.o  energy  barriers  stand  between  that  state  end  steady- state  c  cr~bus 1 1  on. 


Only  considers  t  (.“bust  1  on. 
Sir, liar  to  king's. 


Guide  l< m ration  controlled  by  diffusion  of  o*ycen  through  liquid  oxide  layer. 

Uniform-  particle  t«  *  pera  t ure. 

Vaporwfttion  end  Xl  *•  f  f  us  ion  considered  as  srrips  resistances. 

HpG  r*";;,es  f^C?3  by  diffusion  Jimtcd  reaction  (endrt brnr.ic). 

Heal  gam/lcss  by  conviction,  radiation. 

Tartial  pressure  of  H?0p  negligible  in  the  free  str(am. 

Particle  is  stationary  with  respect  to  surrounding  Casses. 

Paper  only  deals  with  ignition. 

Initial  £^3  layer  thickness  is  uniform. 

Uniform  perticlt  1 i-~  pc  r  a  lure  . 

Inc  icnvective  heat  transfer  coc  f  f  ic  ient ,  h,  is  known. 

Surroundings  treated  as  black  at  temperature  TR;n. 

Prior  to  c*ide  removal,  boron  reaction  rate  is  diffusion  controlled  and  ferr-.s 
P7O3  (»?{,/ ction  follrws  ar.  Arrhenius  type  relation. 

Particle  1 1  •  ;  » r« 1 ure  rcairs  constant  dur  i  ng  r..el  t  i  ng  and  boiling. 

Variation  in  h,  surface  ruction  rate  and  oxide  evaporation  rate  are  r. (elected  during  oxide 
toi ling. 

Water  sapor  rid. kps  the  activation  energy  of  oxide  rvc;  r*r  a  t  ion. 

Initial  particle  ti'perature  is  below  oxide  boiling  ar.d  t  c*r  on  .-.citing. 
f-pOj  fv'-'ed  during  boiling  r.ay  be  neglected. 

So  r.e  t  rad’.ativf  ♦•che'-gf  between  particle  and  cloud. 

Convective  brat  transfer  reelected  in  final  computations. 

The  particle  has  been  heated  to  r  q ui  1  ibr  i irr.  tf-.perature. 
reaction  occurs  at  the  surface. 

The  diffusion  coefficients  for  the  oxidant  and  reaction  products  ere  identical. 

►adiaticn  can  be  neglected  to  a  first  <■;  ~.r  c  >  ir.a  t  Ion  (d  ;  ICOym). 

Particle  te-pi rature  is  above  oxide  boiling  and  boron  rfltlng  prints. 


G 1  f  f  ls  i  r  ri  c  ( •-  trolled  E  *♦  3/4  0?  -*  1/?  F2O3. 
liqjtd  layer  Is  i  sot h,  r— .a  1 . 

1  df  a  1  c  as  . 

Jhsc'ptivf  ♦  q  „  1 1  i  t-r  i  u-  for  Op  at  C  and  IP  steady  stage  dtffu'Wn  cf  0-  within  liquid  layer. 
Plan*  f  -  d  1  “  .  i  4  c  r.al  sl«  b  r-.-drl. 

fio  o»*de  layer. 

Uniform-  particle  tc;f nature. 

Surface  (r'buslirn  i  •■ratios. 
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NOMENCLATURE  FOR  EXPRESSIONS  SUMMARIZED 
IN  TABLE  3 


5.1  KING'S 

CPB1 

CPBS 

S°2 

f 

h 

K 

(mw)b 

(mh)b2o2 

Nu 

P 

P°2 

surf 

po 

B2°3 

PH.,0,  Surf 
PHB02,  surf 

pe2o3,  - 

°RX 

QRX2 

R 

R„ 


RAD 

TS 

V. 

J 

X 


e?03,N2 

DhB02,N2 

V'n2 


AH 


*hvap 


MODEL 

liquid  boron  heat  capacity  (ca1/gm°K) 
solid  boron  heat  capacity  (cal/gm‘>K) 
liquid  boron  oxide  heat  capacity  (cal/gm°K) 

fraction  of  boron  in  the  liquid  phase 

2 

gas-particle  heat  transfer  coefficient  (cal/cm  sec^K) 

mass  transfer  coefficient  for  transport  of  boric  oxide  gas  from 

particle  to  free-stream  (gm-mol/cm  atm  sec) 

boron  atomic  weight  (gm/gm-mol) 

boric  oxide  molecular  weight  (grr./gm-mol ) 

Nusselt  Humber 
total  pressure  (atm) 

oxygen  partial  pressure  in  free  stream  (atm) 

boric  oxide  partial  pressure  adjacent  to  particle  surface  (atm) 
boric  oxide  vapor  pressure  (atm) 

water  gas  partial  pressure  adjacent  to  particle  Surface  (atm) 
HB02partial  pressure  adjacent  to  particle  surface  (atm) 
boric  oxide  partial  pressure  in  free-stream  (atm) 
heat  release  of  B ( s )  +  |o2  -*  2-B203(l)  (cal/gm-mol) 
heat  release  of  B(l)  +  |o2  1 )  (cal/gm-mol) 

gas  law  constant,  82.05  atm-cm^/gm-mol °K 

molar  rate  of  boron  consumption  (gm-mol/sec) 

molar  evaporation  rate  of  boric  oxide  (cm-rr;ol/sec) 

molar  rate  of  removal  of  B2(>3  by  water  reaction  (gm-mol/sec) 

boron  particle  radius  (cm) 

free  stream  gas  temperature  (°K) 

particle  center  temperature  (°K) 

particle  temperature  (°K) 

surroundings  radiation  temperature  (t,K) 

particle  Surface  temperature  (°K) 

molecular  volume  of  species  j  (cm  ) 

oxide  layer  thid.ness  (cm) 

2 

diffusivity  of  gaseous  boric  oxide  in  nitrogen  (cm  /sec) 

2 

diffusivity  of  HB0-,  in  nitrogen  (cm  /sec) 

2 

diffusivity  of  water  gas  in  nitrogen  (err.  /sec) 

heat  absorbed  by  reaction  of  H20  with  ®j03( 1 ) (cal/gm-mol ) 

heat  of  fusion  of  boron  (cal/cm) 

heat  of  vaporization  of  E-,03(  1 )  (cal/gm-mol ) 


TABLE  5.  (Continued) 


Q  evaporation  coefficient  of  boric  oxide  liquid 

Qjj  surroundings  absorptivity 

thermal  diffusivity  of  boron  (cm^/sec) 
o  Stefan-Bol  t2n.ann  constant.  1.35C-101?  cal/cm?sec‘>lc4 

c  particle  emissivl'ty 

v  Herta-Knudsen  impingement  factor  (gm-mol/cm  atm  sec) 

8  time  (sec) 

Pg  boron  density  (gm/cm  ) 

Pg^Q^  boric  oxide  density  (gm/cm  ) 

5.2  MODEL  OF  MEESE  AND  SKIFSTAD 

2 

surface  area,  cm 
specific  heat,  caT/g-°K 
particle  diameter,  cm 
activation  energy  for  evaporation 
heat  of  fusion  of  boron,  cal/g 
heat  of  vapori2ation  of  boron  oxide,  cal/nole 
convective  heat-transfer  coefficient,  cal/cm^-sec-‘>K 
thermal  conductivi ty,  cal/cm- sec-“K 
molecular  weight,  g/g-mole 
mass,  g 
pressure,  atm 
Prandtl  number  of  the  gas 


A 

c 

P 

d 

E 

AHf 

AH 

h 

k 

M 

m 

P 

Pr 


vap 


Q 

R 

RE 

Re 


r 

Ar 

1 

t 


a 

c 


X 

P 

a 


heat  of  combustion  of  boron  with  diatomic  oxygen  to  form  liquid 
cal/niole;  total  energy,  ca? 

molar  rate  of  consumption  of  boron,  o-mole/sec;  gas  constant 
molar  rate  of  evaporation  of  boron  oxide,  gm-mole/sec 
Reynolds  number  based  on  particle  diameter  and  relative  velocity 
between  the  gas  and  particle 
radius,  cm 

oxide  film  thickness,  r^-r^,  cm 
temperature,  °K 
time,  sec 

mole  f paction  of  H^O  present  in  the  gas 

surrounding  absorptivity  (=1.0) 
particle  erissivity  (=0.8) 
particle  burning  rate,  cm/sec 
density,  g/cm 

2  4 

Stefan-Boltrmann  constant,  cal/cm  sec°K 


TABLE  5.  (Continued) 


5.3  MOHAN'S  MODEL 


C 

Ct 

D 

d 

h 

k 

k' 

L 


average  molecular  velocity  normal  to  the  surface  (cm/sec) 

heat  capacity  of  B^O^fL ) (cal/gm"*) 

diffusion  coefficient  (cm^/sec) 

diameter  of  particle  (cm) 

thickness  (cm) 

effective  distribution  coefficient  related  to  the  ratio  of  the 
mole  fraction  of  0^  absorbed  in  B^Ojft)  to  gas  phase  mole  fraction 
distribution  ratio  (liquid  phase  to  gas  phase  0?  concentration) 

sum  of  heat  losses  by  radiation  and  conduction  — -!j— — ■ 

cmc  sec 


P 

qB 

qC 

RB  or  C 
R° 

I 

t 

W. 

1 

X 


total  ambient  pressure  (atm) 

heat  liberated  per  unit  mass  of  boron  consumed  at  B  (cal/g) 

heat  absorbed  per  unit  mass  of  B^O^  vaporized  at  C  (cal/g) 

regression  rate  of  the  surface  (cm/sec) 

universal  gas  constant  (8.316-10  ^erg/gm-mole  "K) 

temperature  (6K) 

time  (sec) 

atomic  weight 

mole  fraction  of  0^  in  the  ambient  atmosphere 

ratio  equilibrium  vapor  pressure  of  B^O^  at  Tc  to  total  pressure 


o  evaporation  coefficient 

c  emissivity 

X  gas  phase  thermal  conductivity  (cal/cm-sec  “K) 

■> 

PB>C  density  boron  and  B^  (gm/cm  ) 

t  time  for  the  new  combustion  regime  to  appear  (sec) 

B.C  surfaces 


5. A 


D 

ddo 

Hi 

iHi 


P. 

i 


VOVCHUK'S  MODEL 

molar  specific  heat 

diffusion  coefficient 

initial  diameter  of  the  boron  droplet 

total  enthalpy  of  species  i 

heat  of  reaction  of  species  i 

molar  flux  of  species  i 

mass  flow  of  the  boron  reacting  in  unit  time 
total  pressure 

partial  pressure  of  species  i 

radiative  heat  flux 

universal  gas  constant 

radius 

temperature 

emissi vi ty 
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x 


o 


thermal  cor, due ti vi ty  of  gas 
molecular  weight 
stoichiometric  coefficients 
density 

Stefan’s  constant 
droplet  combustion  time 


Subscripts 

8  (molten)  boron 

K  particle  or  particle  surface 

•  infinity 

Note:  units  not  stated,  therefore  ary  consistent  set  of  units  may  be  used. 


5.5  EDELMAN’S  MODEL 


N2 


prod 


V'B(Tp) 


B.S 

V 

V 


number  of  moles  of  oxygen  per  mole  of  product 
the  total  (dimensionless)  particle  consumption  rate 
number  of  moles  of  boron  per  mole  of  product 
drag  coefficient 

specific  heat  of  the  gas * (f t?/secZ  °K) 

2  2 

specific  heat  of  boron  at  the  surface  (ft  /sec  °K) 

heat  of  reaction  (calorie/gram) 
surface  reaction  rate  constant 
latent  heat  of  vaporization 
molecular  weight  of  boron  (g/g-mole) 
molecular  weight  of  nitrogen  (g/g-mole) 
molecular  weight  of  oxygen  (g/g-mole) 

molecular  weight  of  the  product  (g/g-mole) 
average  molecular  weight  at  the  particle  surface  (g/g-mole) 
pressure  (in  atmospheres) 
vapor  pressure  of  boron  (atm) 

particle  radius  (feet) 
gas  temperature  (°K) 
particle  temperature  ("K) 
temperature  of  the  wall  ( ° K) 
gas  velocity  (ft/sec) 
particle  velocity  (ft/sec) 

mass  fraction  of  boron  at  the  surface  of  the  particle 
mass  fraction  of  nitrogen  at  the  surface 

mass  fraction  of  oxygen  at  t fie  surface  of  the  [article 
rass  fi action  of  oxygen  in  the  Surrounding  (gas) 


TABLE  5.  (Continued) 


Yprod,S 

a 


c 

X 


mass  fraction  of  the  product  at  the  particle  surface 
fraction  of  vaporized  product 
particle  bulk  density  (slugs/ft^) 
emi ssi vi ty 

the  thermal  conductivity  of  nitrogen  (lbs  force/sec  °K) 
gas  density  (slugs/ft^) 

density  at  the  surface  of  the  particle  (slugs/ft  ) 
Stefan  Boltzmann  constant  (lbs  force/ft-sec  °K^) 


o 


and  mass  transfer.  In  most  of  the  formulations  the  energy  equation  is  the 
central  element  of  the  model,  with  the  size  history  and  other  transport  re¬ 
lationships  representing  auxiliary  elements. 

King  (Refs.  11,  12)  considers  three  stages  in  the  boron  combustion 
process:  below  the  phase  change  temperature  of  the  boron,  at  the  point  of  boron 
phase  change,  and  above  the  temperature  at  which  phase  change  occurs  (2450°K). 
Meese  and  Skifstad  (Ref.  13),  who  based  their  model  on  King's  1972  paper 
(Ref.  12)  also  include  boil-off  of  the  oxide  layer.  In  King's  model  it  is 
assumed  that  the  particle  is  moving  at  the  velocity  of  the  surrounding  gas, 
while  Meese  and  Skifstad  and  Edelman  et  al .  allow  for  slip  between  particle 
and  gas  phase  through  the  use  of  a  convective  heat  transfer  coefficient. 
Vovchuck  (Ref.  14)  does  not  include  convective  heat  transfer,  so  that  in 
Table  3  this  section  is  omitted  in  summarizing  Vovchuck's  work,  (A  similar 
approach  is  followed  in  each  case  in  which  a  given  section  of  Table  3  does 
not  apply  to  a  particular  model.)  Particle-gas  radiation  is  included  in 
each  of  the  models  studied,  and  the  expressions  used  in  each  model  for 
this  component  of  the  overall  heat  transfer  process  are  essentially  inter¬ 
changeable,  with  only  those  differences  related  to  specific  model  geometry. 

It  should  be  noted  that  the  radiation  term  shown  for  Vovchuk’s  model  is 
eventually  neglected. 

King  (Refs.  11  and  12),  Meese  and  Skifstad  (Ref.  13),  and  Mohan 
and  Williams  (Ref.  15)  all  consider  oxide  formation  via  the  global  reaction 
B  +  3/4  0^ 1/2  Vovchuk  (Ref.  14)  allows  for  the  formation  of  the 

lower  oxides  BO,  BO2,  and  B202  at  the  particle  surface.  Edelman  el  al . 

(Refs.  16,  17)  considers  the  product  species  to  be  either  BO  or  B2O3, 
but  if  BO  forms  it  does  not  oxidize  further.  Except  for  the  work  of 
Edelman,  et  al . ,  all  of  the  models  assume  that  the  reaction  rate  is  diffu¬ 
sion  controlled.  The  model  described  by  Edelman,  et  al  .  involves  a  finite- 
rate  reaction  of  B  to  either  BO  or  B2O3  (vapor).  King,  and  Meese  and 
Skifstad,  both  include  oxide  consumption  in  their  models,  through  the 
effect  of  water  on  boric  oxide  ( B ^0^  +  E^O  -*■  2  HBO2).  The  only  other 
mechanism  considered  in  any  model  for  oxide  removal  is  evaporation. 
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The  diffusion  equations  used  in  King's  work  are  algebraic  in  form 
and  depend  upon  the  absolute  value  of  the  pressure  at  the  surface  and  at  in¬ 
finity.  Vovchuk  integrates  the  partial  pressure  from  the  particle  surface 
to  infinity.  For  diffusion  through  the  oxide  layer  King  and  Meese  and 
Skifstad,  use  similar  equations  which  vary  only  in  one  constant;  the  con¬ 
stant  used  by  Meese  and  Skifstad  was  derived  from  data  unavailable  to 
King  and  thus  provides  a  better  fit  of  the  available  experimental  data. 

A  major  feature  of  King's  model  (Refs.  11,  12)  is  the  expression 
representing  boron  phase  change;  Meese  and  Skifstad  (Ref.  13)  use  a  similar 
equation,  derived,  however,  in  a  slightly  different  fashion.  Edelman,  el  al . 
(Refs.  16,  17)  provide  an  expression  for  the  finite-rate  evaporation  of  boron. 
Phase  change  does  not  directly  enter  the  models  reported  by  Vovchuk  (Ref.  14) 
or  Mohan  and  Williams  (Ref.  15).  King,  and  Meese  and  Skifstad,  again  use  a 
similar  expression  for  boric  oxide  evaporation,  while  Vovchuk,  and  Edelman, 
et  al .  have  no  oxide  layer  to  remove.  The  oxide  phase  change  model  incorpor¬ 
ated  by  Mohan  and  Williams  instead  of  being  equilibrium  controlled,  is 
finite  rate.  As  noted  earlier,  Edelman,  et  al .  describe  the  only  model  which 
includes  finite-rate  kinetics,  but  in  this  model  it  is  assumed  that  the  par¬ 
ticle  has  ignited  and  no  B203  layer  is  present,  so  that  oxidation  of  the 
boron  proceeds  directly  to  either  BO  or  Vovchuk  considers  only  boron 

combustion  but  in  the  diffusion  controlled  mode. 

The  end  product  of  Vovchuk's  model  is  an  expression  for  the  particle 
burning  time,  and  an  expression  for  burning  time  is  also  obtained  by  Mohan 
and  Williams.  In  the  model  described  by  Meese  and  Skifstad,  it  is  assumed 
that  no  other  phenomena  occur  while  the  oxide  layer  boils  off,  leading  to  an 
independent  expression  for  the  oxide  layer  boil-off  time.  All  of  the  models 
except  Vovchuk's  provide  expressions  for  the  size  decrease  of  the  particle 
with  time  as  combustion  occurs;  for  those  models  which  incorporate  an  oxide 
layer  expressions  for  the  time-rate  of  change  of  the  oxide  layer  thickness 
are  also  defined  by  the  analysis.  In  three  of  the  models:  those  of  King, 
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Meese  and  Skifstad,  and  Mohan  and  Williams,  ignition  time  is  also  provided. 

In  each  case  ignition  is  defined  as  the  point  at  which  the  oxide  layer  boils 
or  flashes  off. 

What  this  review  indicates  is  that  the  two  most  potentially  useful 
models  are  those  described  by  King  (Refs.  11,  12)  and  by  Edelman,  et  al . 

(Refs.  16,  17).  Each  of  these  approaches  is  the  most  general  of  its  respec¬ 
tive  formulation;  diffusion  controlled  bur-ing  in  the  King  approach,  in¬ 
cluding  the  effects  of  an  oxide  layer,  and  combined  diffusion  and  finite- 
rate  kinetics,  but  without  an  oxide  layer,  on  the  model  of  Edelman,  et  al . 

The  other  approaches  are  either  a  derivative  of  these  more  general  models  or 
considerably  more  restricted  in  application. 

Because  it  includes  both  diffusion  and  finite-rate  chemical  kinet¬ 
ics  mechanisms,  the  model  proposed  by  Edelman,  et  al.  is  useful  for  examining 
the  available  data  for  boron  particle  consumption  rate  as  a  function  of  part¬ 
icle  size.  Since  the  oxide  layer  is  assumdd  not  to  exist,  this  approach  can 
be  interpreted  as  a  model  for  the  combustion  process  after  particle  ignition. 
It  is  thus  useful  for  examining  kinetics  effects  in  boron  consumption,  keep¬ 
ing  in  mind  that  the  results  obtained  represent  a  limiting  case  of  the  overall 
boron  combustion  process. 

2.2.2  Mo  del  Descrip t ion 

The  boron  combustion  model  considered  in  this  work  involves  the  sim¬ 
ultaneous  solution  of  these  coupled  differential  equations  describing  the  vel¬ 
ocity,  size,  and  temperature  of  the  particle.  These  equations  can  be  written: 


where 


Cp  =  drag  coefficient 
pg  =  gas  density  (slugs  per  cubic  foot) 
Vg  =  gas  velocity  (ft/sec) 

Vp  =  particle  velocity  (ft/sec) 

R  =  particle  radius  (feet) 
and  6s  =  particle  bulk  density  (slugs/ft  ) 


dR 

dt 
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where  a  .=  the  total  (dimensionless)  particle  consumption  rate  (see  Eq.  21) 
A  =  the  thermal  conductivity  of  nitrogen  (in  pounds  force/sec  °K) 

Nu  =  Nusselt  number 

2  2 

Cp  =  specific  heat  of  the  gas  (in  ft  /sec  °K),  and  the  remaining 
variables  as  already  defined. 


dTR 
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whe  re 


O  p 

C  =  specific  heat  of  boron  at  the  surface  (ft  /sec  °K) 
Ms 

Tg  =  gas  temperature  (°K) 

Tp  =  particle  temperature  (°K) 

L  =  latent  heat  of  vaporization 

A  =  number  of  moles  of  oxygen  per  mole  of  product 
«s  =  surface  reaction  rate  constant 


3 

=  density  at  the  surface  of  the  particle  (slugs/ft  ) 


Yn  =  mass  fraction  of  oxygen  at  the  surface  of  the  particle 
U2,s 

Mg  =  molecular  weight  of  boron  (g/g-mole) 

aHr  =  heat  of  reaction  (c'alorie/gram) 

Mn  =  molecular  weight  of  oxygen  (g/g-mole) 
u2 


4 

o  =  Stefan-Bol tzmann  Constant  (lb-force/ft-sec  °K  ) 

c  =  emissivity 

Tw  =  temperature  of  the  wall  (°K) 

and  other  parameters  are  as  defined  above. 

Tt.e  values  for  the  remainder  of  the  time-varying  parameters  are  treated  as 
simple  algebraic  equations  (using  the  current  value  of  Tp  and  R)  as  follows 
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where  a  =  fraction  of  vaporized  product 

Pw  D / T  x  =  vapor  pressure  of  boron  (atm) 

’  p ' 

Ms  =  average  molecular  weight  at  the  particle  surface  (g/g-mole) 

Yd  .  =  mass  fraction  of  boron  at  the  surface  of  the  particle 

DO 

P  =  pressure  (in  atmospheres) 

and  other  variables  are  as  previously  defined. 
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where 


a  Yn  mass  fraction  of  oxygen  in  the  surrounding  (gas) 

u  2 » 00 

B  =  number  of  moles  of  boron  per  mole  of  product 

and  the  remaining  variables  are  as  previously  defined. 
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where  ^prod  s  =  msss  fraction  of  the  product  at  the  particle  surface 
Mp^d  =  molecular  weight  of  the  product  g/g-mole) 
and  the  remaining  parameters  are  defined  above 


=  average  molecular  weight  at  the  surface  = 
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where  Mw  =  molecular  weight  of  nitrogen  (g/g-mole) 

Yn  =  mass  fraction  of  nitrogen  at  the  surface 
•  »2  9  ^ 

and  the  remaining  variables  as  previously  defined. 
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In  this  formulation  the  coupling  is  provided  by  the  appearance 
in  the  energy  equation  of  the  parameter  a,  the  total  (dimensionless)  part¬ 
icle  combustion  rate,  which  is  itself  a  function  of  the  surface  concen¬ 
trations  of  both  02  and  B.  Thus,  the  solution  problem  is  one  of  solving 
three  highly  nonlinear,  strongly  coupled  ordinary  differential  equations. 
The  code  developed  for  this  work  implemented  a  Newton-Raphson  scheme  to 
solve  this  coupled  set  of  equations,  and  this  approach  proved  to  be  re¬ 
liable  and  relatively  rapid.  Thus  the  approach  was  then  used  to  investi¬ 
gate  different  aspects  of  the  combustion  of  boron  through  comparison  of 
results  with  available  experimental  data. 

2.2.3  Boron  Combustion  Analyses 

The  model  described  in  this  report  combines  both  diffusive 
transport  and  chemical  kinetic  mechanisms  and  can  thus  be  used  to  examine 
kinetics  limitations  on  boron  combustion  as  a  function  of  particle  size, 
ambient  oxygen  concentrations,  pressure  and  ambient  temperature.  Each  of 
these  variables  is  of  direct  practical  interest  because  each  impacts  the 
use  of  boron  in  a  combustion  process.  For  example,  one  strategy  for  maxi¬ 
mizing  boron  combustion  efficiency  is  to  hold  the  boron  in  an  oxygen-defi¬ 
cient  environment  until  it  is  heated  sufficiently  to  avoid  oxide  layer 
formation.  Under  this  scenario,  the  kinetics  of  boron  consumption  at 
elevated  gas  temperatures  becomes  important.  As  another  example,  in  many 
propulsion  system  applications  a  boron-sl urrv-fuel ed  vehicle  would  be 
required  to  operate  over  a  wide  range  of  altitudes  and  flight  Mach  numbers. 
This  impacts  both  the  temperature  in  the  combustor  and  the  ambient  pres¬ 
sure.  Thus,  these  impacts  on  the  kinetics  of  boron  consumption  are  of 
interest. 

As  particle  size  decreases,  at  constant  pressure,  the  effects  of 

chemical  kinetics  increase  relative  to  those  of  diffusion.  This  is  clearly 

shown  by  the  results  presented  in  Fig.  16.  For  example,  for  1  atm  pressure, 

2 

the  transition  from  kinetics  control  d)  to  diffusion  control  (t^  «  d  ) 

occurs  for  particles  somewhat  larger  than  200  ym.  However,  this  transition 


Figure  16.  Predicted  Ki netic-to-Diffusion 
Controlled  Combustion  of  Boron 


is  itself  a  function  of  pressure,  as  the  resul ts  indicated  in  Fig.  16  also 
show.  At  low  pressure,  typical  of  high-altitude  ramjet  operation,  the 
transition  to  diffusion  control  occurs  at  diameters  so  large  as  to  be  well 
outside  the  range  of  practical  interest  (for  this  hot-gas,  hot  particle  case 
while  at  high  pressures,  such  as  might  be  encountered  in  an  advanced  gas 
turbine  combustor,  the  consumption  process  is  di ffusion-controll ed  even  for 
quite  small  particles. 

In  most  applications,  boron  particle  sizes  must  be  kept  small, 
even  at  elevated  pressures,  to  avoid  long  burn  times.  This  is  illustrated 
by  Fig.  17,  in  which  the  results  depicted  on  Fig.  16  are  cross-plotted  to 
indicate  the  effects  of  pressure  on  consumption  time  f.or  a  given  particle 
size.  It  is  immediately  evident  from  this  figure  that  for  large  particle 
sizes  burn  times  greater  than  10  ms  would  be  expected  (in  the  absence  of 
a  B2O2  coating;  considerably  larger  burn  times  can  be  anticipated  for 
B^O^-clad  particles).  Note  also  that  fop  large  particles  the  burn  time 
tends  to  become  independent  of  pressure  at  high  pressures  as  the  transition 
to  diffusion  control  is  completed.  For  small  particles,  on  the  other  hand, 
kinetics  effects  dominate  throughout  the  pressure  range  studied,  and  at  a 
few  atmospheres  pressure  burn  times  of  a  fraction  of  a  millisecond  can  be 
achieved . 

While  the  data  shown  in  Fig.  16  represent  high-temperature  boron 
combustion,  much  of  the  available  data  involves  ignition  and  combustion  in  a 
cold  environment.  Figs.  18  and  19  provide  data  at  several  pressure  levels, 
for  several  particle  sizes,  at  two  ambient  oxygen  concentrations,  and  show 
comparisons  between  the  data  and  predictions  carried  out  using  the  present 
model.  Overall,  the  agreement  between  the  data  and  the  model  predictions  is 
quite  good,  although  the  model  does  not  represent  well  the  decrease  in  burn 
time  at  pressures  above  about  20  atmospheres  that  may  be  indicated  by  the 
data.  Another  aspect  of  the  model  prediction  is  also  of  interest:  at  low 
pressures,  for  smaller  particles,  the  model  predicts  that  the  boron  consump¬ 
tion  process  will  quench,  leaving  a  non-negl i gi bl e  residue.  This  residue 
particle,  while  not  insignificant  from  a  combustion  inefficiency  standpoint. 


is  not  large  and  could  be  missed  in  experimental  studies  unless  its  presence 
was  specifically  being  investigated.  The  quenching  process  appears  to  occur 
at  high  pressures  for  smaller  particles,  and  could  be  of  considerable  sig¬ 
nificance  in  ramjet  operation,  where  afterburning  can  be  continuing  in  the 
exhaust  nozzle. 

2.2.4  Finite-Rate  BgO^  Condensation 

While  the  existence  of  as  an  oxide  coating  on  a  boron  particle 
is  of  critical  importance  in  the  examination  of  boron  particle  ignition,  the 
condensation  of  the  oxide  out  of  the  gas  phase  is  also  critical  to  the  over¬ 
all  performance  of  a  boron  slurry  fueled  ramjet.  This  follows  from  the  fact 
that  approximately  30%  of  the  available  energy  of  boron  combustion  is  tied 
up  in  the  phase  transformation  between  liquid  and  gaseous  ^2^3'  Where  re- 
condensation  of  B^O^  occurs  --  in  the  combustor  or  in  the  nozzle  --  is  thus 
critical  to  achieving  the  optimum  performance  levels  utilizing  boron  fuels. 

Finite-rate  condensation  effects  can  be  analyzed  using  single  part¬ 
icle  nucleation  and  growth  theory.  For  example,  nucleation  rate  can  be 
approximately  using  the  expression 


J  = 


(27) 


from  classsical  kinetic  theory,  where 

K  is  the  Boltzmann  constant, 

T  is  the  local  temperature,  and 
B  is  a  rate  constant; 

AG  is  the  free  energy  of  formation  of  a  critical  cluster,  itself  a  function 
of  the  surface  tension  of  the  material,  the  work  of  cluster  formation,  and 
the  number  of  molecules  in  a  critical  droplet.  Thus  AG=f-|  (o,  w*,  n*).  w* 
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and  n*  are  both  functions  of  r*,  where  r*  is  the  radius  of  a  critical  drop¬ 
let  A  number  of  expressions  have  been  devised  for  the  terms  representing 
the  work  of  cluster  formation  and  the  number  of  molecules  in  a  critical 
droplet.  For  example,  in  the  expression  developed  by  Yang  (Ref.  18)  and 
used  by  Edelman  and  Economos  (Ref.  19)  to  examine  boron  condensation 
phenomena,  the  work  of  cluster  formation  is  given  by 
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where  o is  the  liquid  surface  tension,  and  the  number  of  molecules  in  a 
critical  drop  is  expressed  as 


n* 
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in  which  is  the  molecular  weight  and  the  density  of  the  liquid 
condensate. 


The  radius  of  a  critical  drop  can  be  expressed  as 
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where  P  is  the  vapor  partial  pressure  and  P  is  the  equi  1  i bri urn -saturation 
partial  pressure,  i.e.,  the  vapor  pressure  of  a  vapor  in  equilibrium  with  a 
plane  surface  of  condensed  phase.  RQ  is  of  course  the  universal  gas  constant. 
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Once  the  nucleation  process  produces  a  critical  size  droplet,  it 
continues  to  grow  by  the  additional  deposition  of  condensing  vapor.  A  growth 
law  which  has  been  used  to  calculate  water  vapor  condensation  phenomena  is 
given  by  the  expression 


,  (P  -  P  )  /  M,  \1/2 

dr  _  v  v  Voo '  I  L  ] 
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Finally,  the  total  production  rate  of  condensed-phase  B^O^  can  be  written 
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where  the  growth  rate  term  is  summed  over  several  classes  of  particles 
defined  by  a  mean  radius  r-,  i.e., 

J 


(33) 


in  which  gives  the  number  density  of  particles  of  radius  class  j. 

Application  of  this  approach  to  the  finite-rate  condensation  of 
requires  knowledge  of  the  surface  tension  of  B^O^  and  its  equilibrium- 
saturation  partial  pressure.  These  data  can  be  estimated  for  the  tempera¬ 
ture  ranges  of  interest  through  modest  extrapolation  of  available  data. 

For  example,  as  shown  by  Fig.  20,  data  included  in  the  JANNAF  tables  allows 
a  modest  extrapolation  of  the  low-temperature  data  to  fit  the  empirical 
express  ion 

Pv„  =  exP  (19-03  -  44f“)  (34) 
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Figure  20.  B~0,  Vapor  Pressure  Data 


Surface  tension  data  requires  a  rather  larger  extrapolation,  as  shown  by  the 
data  included  in  Fig.  21. 

The  physical  data  described  here  and  the  expressions  for  nucleation 
and  condensation  set  out  in  this  section  were  used  in  the  work  outlined  in 
Ref.  19  to  compute  the  condensation  of  in  a  supersonic  nozzle  expan¬ 
sion  process;  the  condensation  model  was  incorporated  in  an  analysis  which 
coupled  finite-rate  chemistry  for  the  B-O-H-N  system  with  finite-rate  con¬ 
densation  phenomena.  Fig.  22  shows  the  nozzle  pressure  distribution  data 
obtained  in  a  test  involving  condensation.  These  experimental  data 
show  a  definite  departure  from  the  otherwise  isentropic  expansion  associated 
with  the  conical  nozzle  geometry.  Also  given  in  this  figure  are  the  results 
of  a  calculation  based  on  the  assumption  that  the  latent  heat  of  con¬ 
densation  is  released  i nstantaneously ,  resulting  in  a  Rayleigh  line  process 
assumed  to  occur  in  the  vicinity  of  the  condensation-induced  pressure  rise 

as  indicated  by  the  data.  The  heat  release  used  in  the  Rayleigh  model  is 

» 

based  on  the  condensation  of  all  of  the  available  in  the  flow.  These 
results  indicate  that  the  agreement  of  the  Rayleigh  analysis  and  the  ex¬ 
perimental  data  is  in  fact  quite  good,  although  the  ideal  model  somewhat 
overestimates  the  magnitude  of  the  experimental  pressure  rise.  Cn  the 
other  hand,  the  coupled  finite-rate  condensation  and  gas  phase  kinetics 
approach  agrees  fairly  well  with  the  data  though  the  initial  pressure  change 
is  somewhat  underpredicted .  This  initial  underprediction  casts  some  doubt 
on  the  nucleation  and  growth  laws  used  in  this  study,  but  given  the  meager 
data  available  both  for  the  physical  parameters  which  enter  the  models  and 
for  testing  and  development  of  the  expressions  themselves,  such  a  deviation 
from  measured  data  as  shown  here  is  not  unexpected.  This  work  does,  on  the 
other  hand,  provide  evidence  that  the  effects  of  B^O.^  condensation  can  be 
included  in  analytical  approaches  intended  to  model  the  performance  of  boron- 
fueled  combustors. 


O  EXPERIMENTAL  DATA: 
NOZZLE  EXPANSION 
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Figure  22.  Effects  of  Condensation  in  a 

Supersonic  Nozzle  Expansion 


The  results  presented  in  this  section  indicate  that  a  coupled 
diffusion  and  finite-rate  kinetics  approach  is  required  to  adequately  examine 
boron  combustion  in  the  context  of  its  use  as  a  fuel  in  ramjet  combustors. 

These  results  show  that  coupling  finite-rate  oxidation  kinetics  with  dif- 
fusional  processes  allows  the  behavior  of  boron  combustion  as  observed  in 
laboratory  experiments  to  be  explained,  in  terms  of  burn  time/particl e  size 
and  pressure  and  temperature  relationships.  An  understanding  of  these  phenom¬ 
ena  is  crucial  to  the  utilization  of  boron  in  combustors.  The  computations 
also  clearly  indicate  quenching  of  the  boron  combustion  process  even  for 
clean  (oxide  free)  particles  for  certain  combinations  of  partical  size,  ambient 
temperature,  and  ambient  oxygen  concentration. 

Just  as  clearly  the  work  described  herein  indicates  the  need  for 
a  number  of  improvements .  For  example,  i  incorporation  in  the  modeling  of 
the  presence  of  an  oxide  layer  on  the  particle  is  needed  to  adequately 
account  for  boron  ignition  phenomena.  The  role  and  potential  importance 
of  suboxide  formation  and  disposition  remains  unclear,  and  more  mech¬ 
anistic  and  detailed  models  of  B-0,  B-O-H,  and  B-O-H-C  kinetics  are  clearly 
required.  Additionally,  the  work  reported  herein  also  shows  that  B2C>3  con¬ 
densation  processes  can  be  modeled,  and  the  impact  of  this  on  the  overall 
heat  release  obtained  with  boron  combustion  is  such  that  condensation 

needs  to  be  controlled  to  optimize  combustion  efficiency.  Further  work 
remains  to  establish  appropriate  analyses  of  these  processes.  Finally, 
cloud  and  slurry  spray  combustion  analyses  need  to  be  established  to  aid 
in  the  practical  use  of  boron  slurry  fuels.  The  individual  parts  of  such 
an  analysis  exist,  as  indicated  by  the  spray  modeling  work  considered  in 
the  preceding  section,  the  single-particle  boron  modeling  considered  here, 
and  the  modular  model  work  to  be  discussed  in  the  next  section.  What  is 
required  is  the  merging  of  the  efforts  in  these  three  research  areas. 


2.3 


THE  MODULAR  MODEL  FOR  RAMJET  PERFORMANCE  ASSESSMENT 


The  modular  model  for  ramjet  performance  assessment  represents 
the  most  advanced  and  detailed  combustor  model  applicable  to  sudden- 
expansion  ramjets.  It  has  been  developed  because  the  complex  flame  stabil¬ 
ization,  flame  propagation,  and  spray  combustion  problems  which  have  been 
encountered  in  the  development  of  dump  combustors  require  a  detailed 
analytical  capability  for  their  resolution.  The  ability  to  compute  in  some 
detail  combustion  chamber  flowfields  is  necessary  in  order  to  understand 
the  phenomena  that  occur  in  existing  combustors  and  to  predict  the  perfor¬ 
mance  of  new  combustor  concepts. 

In  its  basic  form,  the  modular  model  approach  has  been  described 
in  a  variety  of  reports  and  papers,  of  which  Refs.  20  and  21  are  examples. 

The  model  described  in  this  report,  however,  differs  from  thore  previously 
described  in  two  key  areas:  first,  the  well-stirred  reactor  cnemistry 
model  is  fully  integrated  with  the  chemistry  modeling  applied  to  i h e  directed 
flow,  and  second,  the  advanced  quasiglobal  model  for  hydrocarbon  combustion 
kinetics  is  incorporated  in  the  code.  This  latter  difference  between  the 
present  model  and  prior  versions  of  the  modular  approach  is  ext  •  mely 
important  because  it  is  only  with  the  advanced  quasiglobal  model  that  it 
is  possible  to  accurately  characterize  the  combustion  kinetics  and  heat 
release  processes  that  occur  in  fuel-rich  regions.  Since  rucl-rich  regions 
make  up  a  significant  portion  of  nonpremixed  ramjet  combustor  flowfields, 
and  can  be  a  significant  cause  of  combustion  inefficiency  in  ramjet  opera¬ 
tion,  this  fuel-rich  combustion  characterization  ability  s  a  key  require¬ 
ment  for  an  overall  ramjet  combustor  pei  fc  rmance  mod  1. 

The  modular  concept  is  that  the  combustor  flowfield  can  be  broken 
down  into  three  major  components,  as  reprt  e>‘ed  schematically  o  rig.  23. 
These  components  are  a  directed  flow  (which  is  computationally  modeled 
as  a  flowfield  satisfying  the  boundary  layer  approximation  that  axial 
gradients  are  negligible  compared  to  radial  gradients,  and  the  static  pres¬ 
sure  is  radially  uniform);  a  reci rculation  zone  (assumed  to  be  represent¬ 
able  as  a  well-stirred  reactor  or  a  combination  of  well-stirred  reactors); 
and  a  turbulent  shear  layer  region,  located  along  the  dividinc  streamline 
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which  separates  the  other  two  regions.  The  shear  layer  serves  to  couple  the 
other  two  model  components;  fluxes  of  species  and  energy  across  this  shear 
layer  form  the  boundary  conditions  on  the  two  computational  regions.  The 
directed  flow  is  assumed  to  be  fully  turbulent,  with  the  turbulent  effective 
viscosity  defined  by  a  conventional  two-equation  turbulence  model  (Ref.  22). 

2.3.1  Stirred  Reactor  Modeling  of  Dump  Combustor  Reci rcul a tion  Zones 


The  well-stirred  reactor,  shown  schematical ly  in  Fig.  24,  is  a 
laboratory  device  in  which  very  high  mixing  rates  are  achieved.  In  general, 
laboratory  stirred  reactors  are  designed  to  ensure  that  the  mixture  within  the 
reactor  is  spatially  uniform,  so  that  unreacted  feed  material  is  continuously 
and  uniformly  mixed  with  combustion  products,  reacting  for  a  time  defined  by 
the  average  residence  time  of  the  reactor  before  exiting.  For  the  limit  of 
perfect  stirring,  this  average  residence  time  is  given  simply  by  the  ratio 
of  the  product  of  the  stirred  reactor  yolyme  and  density  to  the  feed  mass 
flow,  pV/H;  the  product  mass  flux  is  of  course  equal  to  the  feed  mass  flu:. 

In  this  limit  the  governing  equations  for  the  stirred  reactor  state  reduce 
to  algebraic  relations,  which  allows  the  use  of  rapid  and  efficient  solution 
procedures . 

In  the  lower  half  of  Fig.  24,  a  typical  reci.rcul ation  region  in  a 
sudden-expansion  combustor  is  sketched.  Because  of  the  high  turbulence 
intensity  and  large  mixing  rates  generated  in  the  reci rcul  ation  zone,  the 
limiting  behavior  of  this  region  can  be  thought  of  as  well-mixed,  that  is, 
the  state  of  the  recirculation  region  is  defined  by  the  chemical  kinetics 
of  the  reactions  occurring  and  not  by  the  mixing  rate.  In  this  sense  the 
volume  within  the  reci rculation  region  sketched  in  Fig.  24(b)  is  similar  to  the 
volume  within  the  laboratory  stirred  reactor  in  Fig.  24(a),  and  the  same 
solution  technique  can  be  used  to  obtain  the  thermo  chemical  static.  A  com¬ 
parison  of  Figs.  24(a)  and  24(b)  also  show's  the  major  difference  between  the 
recirculation  region  and  the  laboratory  stirred  reactor.  In  the  laboratory 
stirred  reactor,  discrete  reactant  and  product  streams  can  be  identified,  but 
in  the  well-stirred  reactor  model  of  a  sudden-expansion  recirculation  zone, 
separate  reactant  and  product  streams  cannot  be  defined.  Instead,  reactants  enter 
and  pioducts  leave  the  recirculation  region  by  turbulent  diffusion  through 


the  shear  layer  which  surrounds  the  dividing  streamline.  Thus,  the  re¬ 
actant  stream  is  the  flux  of  reactants  passing  through  the  shear  layer, 
integrated  over  the  surface  area  of  the  dividing  streamline,  and  the  pro¬ 
duct  stream  is  defined  similarly.  These  fluxes  are  of  course  equal  and 
oppos i te. 

The  equations  describing  the  transport  of  energy  and  species  in 
the  well-stirred  reactor  reduce  to  the  following  relations: 


Conti nui ty 


•I  -o 

m  =  m 


Spec i es , 


=  m.  +  Vw. 
dt  i  i 


Energy, 


h  =  ?  h!ai  =  hI  +  Q/m 


Note  that  in  this  set  of  equations  the  species  transport  equation  is 
written  in  nonsteady  form.  This  formulation  has  been  adopted  to  facilitate 
solution  of  the  stirred  reactor  governing  equations  with  finite-rate  chemical 
kinetics;  the  steady-state  stirred  reactor  solution  is  obtained  when 
da.j/dt  -*■  0. 

Equations  (35-37),  along  with  expressions  for  the  volumetric  species 
production  rates,  for  the  enthalpy  of  individual  species  as  a  function  of 
temperature,  and  the  equation  of  state,  define  the  temperature  and  species 
concentrations  in  the  stirred  reactor,  given  the  inflow  rates  for  species 
and  enthalpy.  In  the  modular  model,  the  net  inflow  of  species  and  enthalpy 
are  each  expressed  as  line  integrals  involving  gradients  evaluated  along  the 
dividing  streamline,  so  that  the  energy  and  species  conservation  equations 
for  the  stirred  reactor  are  written,  respectively. 


and 


where  pvy  is  evaluated  using  the  outer  flowfield  solution  in  the  region  of  the 
dividing  streamline. 

Because  in  the  modular  mode!  me  teed  rates  into  the  recirculation 
region  are  defined  by  the  fluxes  of  reactants  and  products  through  the 
shear  layer,  this  region  becomes  a  key  element  of  the  model.  Nevertheless, 
the  shear  layer  is  modeled  simply  as  a  region  of  linear  gradients;  that  is, 
the  gradient  in  a  quantity  <P  at  the  dividing  streamline  is  approximated  by 


Moreover,  the  width  of  the  shear  layer  is  assumed  to  be  given  by  the  linear 
expression 


1  =  ax  +  b  (41 ) 

in  which  a  and  b  are  defined  through  comparison  of  model  predictions  with 
experimental  data.  The  results  obtained  with  the  modular  approach  have 
been  found  to  be  relatively  insensitive  to  the  values  of  the  coefficients 
a  and  b  in  Eq.  (41);  for  the  results  presented  in  this  report,  a  =  0.03 
and  b  =  0.  07 . 

2.3.2  Parabol ic  Mixing:  The  Director-Flow  Model 

The  second  major  element  of  the  modular  model  for  a  sudden- 
expansion  combustor  is  the  formulation  for  the  di rected-fl ow  portion  of  the 
combustor  flowfield.  It  is  assumed  that  the  boundary-1 ayer  approximations 
apply  to  this  part  of  the  flowfield.  For  a  steady,  axisymmetric  flow,  the 
describing  equations  may  be  written  as  follows: 

Global  Continuity, 


9r_pu  3rpv 

3x  3r 


Species  Diffusion  for  the  ith  Specie, 


(42) 


pu 


+  pv 


Momentum  Equation, 


(43) 
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3u 

3x 
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Energy  Equation, 


3H  3H  1  3 

pU  3^  +  pV  3?  =  r  37 


rpv. 


Pr 


a»  - 1  p-n  -  i) 

3r  \ Sc  j 


-  ”  £  (t) 

These  equations,  along  with  expressions  for  the  enthalpy 


H  =  h  +  (u^/2)  and  h  =  (T) 

i 

and  the  equation  of  state 


pRT  £  (a‘/W,) 


r  v 


(45 


(46 


(47 


can  be  solved,  given  an  expression  for  the  turbulent  eddy  viscosity  py  =  Pvy 
This  is  defined  by  a  two-equation  turbulent  kinetic  energy  model,  which 
introduces  transport  equations  for  the  turbulent  kinetic  energy  and  its 
dissipation  rate.  In  boundary-layer  form,  these  equations  can  be  written 
as  follows: 


Turbulent  Kinetic  Energy 
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Turbulence  Energy  Dissipation, 
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where  Uy  =  C^pk^/e.  A  standard  set  of  two-equation  model  coefficients  has  been 
used  for  all  of  the  computations  described  in  this  report:  thus  Cy  =  0.09. 
C£1  =  1.40,  C£2  -  1.95,  ok  =  1.00,  and  =  1.22. 

2.3.3  Chemical  Kinetics:  The  Advanced  Quasiglobal  Model 

It  has  long  been  noted  that  an  important  limitation  of  the  origi¬ 
nal  basic  quasiglobal  formulation  is  its  inapplicability  under  fuel-rich 
conditions;  e  g.,  for  fuel-air  ratios  greater  than  stoichiometric.  For 
equivalence  ratios  greater  than  about  1.  5,  the  adiabatic  flame  temperature 
predicted  by  the  basic  quasiglobal  approach  is  greater  than  that  observed 
in  experiments,  and  the  discrepancy  between  the  predicted  and  measured 
adiabatic  flame  temperatures  becomes  larger  as  the  combustion  process  be¬ 
comes  more  fuel-rich.  This  is  a  serious  difficulty  with  respect  to  diffu¬ 
sion  flames  in  general  and  in  liquid-fueled  sudden-expansion  combustor  ram¬ 
jet  applications  in  particular.  In  these  devices  the  local  fuel -air  ratio 
in  and  near  the  flame-stabilizing  recirculation  regions  has  been  found  to 
be  well  above  stoichiometric.  Gas  turbines  as  well  often  have  fuel-rich 
primary  zones  in  their  combustion  chambers.  Thus  a  clear  need  exists  for 
an  improved  quasiglobal  model,  in  order  to  examine  kinetic  processes  in 
fuel-rich  combustor  regions. 


The  basic  reason  that  the  original  quasiglobal  model  formula¬ 
tion  fails  in  the  prediction  of  combustion  processes  in  fuel-rich  regions  is 
that  it  does  not  include  the  fuel  pyrolysis  steps  that  occur  in  high  temp¬ 
erature,  oxygen  deficient  regions,  but  instead  assumes  a  pure  oxidation 
step  for  fuel  consumption.  In  the  fuel  rich  regions,  both  pure  pyrolysis 
and  oxidative  pyrolysis  steps  occur  to  lead  to  the  breakdown  of  the  fuel 
into  lower  molecular  weight  components.  At  the  same  time,  soot  formation 
occurs:  while  this  process  does  not  soak  up  a  large  amount  of  energy,  it 
does  result  in  the  formation  of  particulates  that  strongly  increase  the 
radiative  heat  loss  from  the  flame.  This  radiation  enhancement  results  in 
both  a  decrease  in  the  overall  temperature  (and  thus  can  directly  affect 
the  rate  of  heat  release)  and  an  increase  in  the  heat  load  to  the  combustor 
walls.  Eased  on  these  obserya tions ,  the  advanced  quasiglobal  model 
incorporates  both  fuel  pyrolysis  and  soot  formation  steps. 
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In  the  advanced  quasiglobal  kinetics  approach,  the  fuel  to 
be  modeled  is  broken  down  into  aromatic  and  aliphatic  components:  soot 
formation  is  dependent  on  the  amount  of  aromatic  fuel  present.  The 
aromatic  and  aliphatic  components  are  modeled  using  rates  which  have 
been  developed  for  toluene  and  iso-octane,  respectively.  That  is,  the 
model  fuel  component  which  can  be  treated  as  aromatic  is  assigned  rates 
which  have  been  developed  by  extensive  comparison  of  toluene  predictions 
with  experimental  data,  and  the  model  fuel  component  that  can  be  treated 
as  aliphatic  is  assigned  the  rates  developed  for  iso-octane.  Each  of 
these  fuels  break  down  into  secondary  fuels,  C2H2  from  toluene  and  C2H4 
from  iso-octane,  following  the  pure  pyrolysis  pathway,  as  well  as  par¬ 
tially  oxidizing  by  attack  from  both  0 ^  an^  OH.  The  secondary  fuels  then 
react  with  the  remaining  OH  and  O2,  and  the  wet  CO  mechanism  follows. 

Table  6  illustrates  these  reactions,  along  with  soot  formation/consump- 
tion  reaction. 

While  this  remains  a  formidable  list  of  reactions,  our  studies 
show  .hat  the  accuracy  of  the  predictions  is  most  sensitive  to  the  wet-CO 
portion  of  the  mechanism.  This  observation  indicates  that  further  simpli¬ 
fications  are  possible  depending  on  the  information  being  sought.  For 
example.  Table  7  shows  a  simplified  mechansim  intended  to  examine  the  heat 
release  distribution  in  a  ramjet  combustor  where  the  flame  is  stabilized  by 
a  recirculation  zone.  In  this  example,  many  of  the  intermediates  important 
in  a  lower  temperature  ignition  process  are  presumed  to  be  negligible  because 
of  the  presence  of  the  recirculation  zone.  Also,  in  this  example  the  fuel, 
JP-7,  has  been  modeled  as  £-^24*  anc*  ^as  ^een  treated  as  being  aliphatic  in 
character. 

2.3.4  Fuel  Jnject ion  Model 

In  addition  to  the  basic  components  of  a  parabolic,  directed  flow 
analysis  and  a  well-stirred  reactor  formulation,  the  modular  concept  in¬ 
cludes  other  elements  of  the  dump  combustor  flowfield,  for  example,  the 
fuel  injection  process.  The  detail  of  the  computation  provided  by  the  use 
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of  a  parabolic  directed  flow  analysis  is  the  key  feature  of  the  modular 
model  that  allows  the  inclusion  of  a  fuel  injection  module.  This  is 
particularly  true  in  the  case  of  liquid  fuel  injection,  for  at  the  fuel/air 
ratios  appropriate  for  ramjet  combustor  operation,  the  liquid  fuel  streams 
initially  occupy  a  very  small  portion  of  the  overall  combustor  cross- 
sectional  area.  The  liquid  fuel  injection  model  makes  use  of  a  combina¬ 
tion  of  empirical  information  and  turbul ent  mixing  calculations.  For  ex¬ 
ample,  the  fuel  jet  penetration  from  the  wall  is  computed  through  the  use 
of  an  empirical  penetration  correlation  (Ref.  23),  using  a  breakup  time 
correlation  (Ref.  24)  to  estimate  the  downstream  distance  at  which  penetra- 
.tion  is  to  be  computed.  That  is,  it  is  assumed  that  the  fuel  jet  has  turned 
parallel  to  the  airflow  at  the  axial  position  at  which  the  initial  fuel  jet 
has  broken  up  into  droplets,  as  given  by  the  breakup  time  correlation  and 
the  local  airflow  velocity. 

Since  the  basic  modular  model  formulation  involves  an  axisymmetric 
flowfield,  individual  fuel  jets  cannot  be  resolved,  and  it  is  assumed  that 
the  liquid  fuel  spray  forms  an  annulus  whose  cross-sectional  area  may  either 
be  specified  or  computed  based  on  an  assumed  fuel  spray  bulk  velocity.  A 
bulk  spray  evaporation  correlation  is  then  used  to  compute  the  fuel  vapor¬ 
ization  rate;  this  correlation  (Ref.  25)  is  a  function  of  the  initial  velo¬ 
city  and  temperature  difference  between  the  fuel  spray  and  the  surrounding 
air  stream.  Spreading  of  the  fuel  jet  is  computed  through  use  of  a  turbu¬ 
lent  mixing  hypothesis  as  for  the  mixing  process  in  the  remainder  of  the 
parabolic  flow.  Fig.  25  shows  the  results  of  a  computation  of  the  fuel 
injection  process  for  three  fuel  injectors,  located  in  the  combustor  inlet 
wall,  along  with  the  centerline,  and  in  a  midstream  position.  Shown  are 
the  computed  contours  of  the  fuel  mass  fraction,  Op,  with  the  vapor-phase 
fuel  shown  as  the  solid  line  and  the  liquid-phase  fuel  as  the  dotted  line, 
as  a  function  of  both  axial  and  radial  position  in  the  combustor  inlet. 

For  these  calculations  the  air  inlet  velocity  was  approximately  700  ft/sec 
at  a  temperature  of  1600°K;  the  overall  fuel/air  equivalence  ratio  was 
0.6.  The  results  shown  in  Fig.  25  provide  a  good  example  of  the  detail 
of  the  fuel  injection  process  available  through  use  of  tin's  aspect  of 
modular  modeling. 
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2.3.5  It eration  Procedure 

The  overall  logic  of  the  calculation  procedure  is  shown  schemat¬ 
ically  in  Fig.  26.  Inlet  conditions  can  be  defined  upstream  or  downstream 
of  the  fuel  injection  station  as  required;  if  liquid  fuel  injection  is 
specified,  the  penetration  and  spreading  of  the  liquid  fuel  jets  are  com¬ 
puted  as  described  in  the  preceding  section.  The  computation  proceeds  to 
the  dump  plane  in  the  sudden-expansion  combustor,  at  which  point  the  depend¬ 
ent  variables  in  the  flowfield  are  stored  for  later  use  in  the  recirculation 
zone  iteration  procedure. 

An  initial  state  of  the  stirred  reactor  must  be  specified  in  order 
to  begin  the  iteration  procedure.  This  state  is  reasonably  arbitrary,  except 
that  a  reacted  temperature  level  must  be  specified.  The  shape  of  the  divid¬ 
ing  streamline  separating  the  directed,  flow  and  the  recirculation  zone  must 
also  be  specified,  along  with  the  shear  stress  level,  expressed  as  a  "skin 
friction"  coefficient.  The  initial  stirred  reactor  state,  in  conjunction  with 
a  model  for  the  shear  layer  between  the  recirculation  zone  and  the  directed 
flow  then  defines  the  boundary  conditions  for  the  parabolic  directed  flow 
calculation.  Directed  flow  computations  are  carried  out  to  the  axial  station 
at  which  the  end  of  the  reci rcul a ti on  zone  has  been  defined;  as  part  of  these 
computations  the  diffusive  flux  of  species  and  energy  across  the  dividing 
streamline  is  computed.  These  fluxes  then  define  a  new  set  of  stirred  re¬ 
actor  "feed  rates,"  i.e.,  species  and  energy  fluxes  into  the  recirculation 
zone  region.  Convergence  of  the  procedure  is  defined  by  the  change  in  the 
stirred  reactor  feed  rates  from  iteration  to  iteration;  each  species  and 
the  energy  flux  must  change  less  than  1.051  before  the  convercen  e  criteria 
are  satisfied.  If  they  are  not,  the  species  and  energy  fluxes,  and  the  Over¬ 
all  diffusive  mass  flux  and  physical  recirculation  zone  volume  are  used  to 
compute  a  new  stirred  reactor  state.  The  pressure  required  for  this  compu¬ 
tation  is  taken  to  be  the  arithmetic  average  of  the  pressures  computed  as 
part  of  the  directed  flow  solution  at  the  beginning  and  end  of  the  recircu¬ 
lation  region.  The  new  stirred  reactor  state  is  then  used  to  define  the 
directed  flow  boundary  conditions,  and  the  computation  is  restarted  from  the 
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igure  26.  Modular  Model  Calculation  Logic 


dump  station.  When  convergence  is  achieved,  typically  in  3-5  iterations, 
the  combustor  calculation  is  continued  to  the  specified  combustor  exit 
station. 

2.3.6  Modular  Model  Assessment 

A  variety  of  applications  of  different  versions  of  the  modular 
model  have  been  reported  in  the  literature,  Refs.  20  and  21,  and  performance 
calculations  have  been  carried  out  for  several  types  of  sudden  expansion 
ramjet  combustors,  with  inlet  conditions  including  premixed,  uniform  inlet 
profiles,  wall  and  in-stream  orifice  fuel  injectors,  and  wall-mounted 
poppet  type  fuel  injectors. 

An  early  set  of  modular  model  validation  comparisons  included 
the  sudden  expansion  ramjet  combustor  results  described  by  Craig,  et  al. 

(Ref.  26).  This  work  involved  a  combustor  L/D  of  3,  a  dump  area  ratio 
AyA2  of  2.25,  and  an  exit  nozzle  area  ratio  A^/A^  =  0.40.  Cold  flow  and 
reacting  flow  experimental  results  are  both  the  reported  for  this  config¬ 
uration.  In  the  modular  model  computations,  the  cold  flow  data  were  used 
to  evaluate  the  model  coefficient  values  for  the  specific  configuration, 
and  computations  of  a  premixed  dump  combustor  at  a  fuel/air  ratio  of 
0.053  were  then  carried  out  using  these  coefficient  values  for  comparison 
with  the  data  presented  in  Ref.  26.  In  the  experiments,  the  fuel  used 
was  JP-4  which  was  represented  in  the  computations  by  propane.  Both  a 
simple  one-step  global  finite-rate  chemistry  model  and  the  full  hydrocarbon 
oxidation  kinetics  as  represented  by  the  quasiglobal  model  were  used  to 
represent  the  chemical  kinetics  processes  in  the  directed  flow.  A  global 
model  was  used  for  the  recirculation  zone  in  all  cases.  Propane  was  chosen 
to  represent  the  fuel  because  it  adequately  simulates  from  a  heat-release 
standpoint  the  actual  fuel  used  and  because  previous  studies  with  propane 
had  resulted  in  the  development  of  a  one-step  finite-rate  kinetics  model  for 
this  fuel.  This  one-step  model  represents  fairly  accurately  the  ignition 
delay  for  propane-air  (although  not  the  overall  reaction  time)  over  a  range 
of  conditions  of  interest  in  this  work.  Initial  conditions  included  an  inlet 
total  temperature  of  554°K  and  an  inlet  static  pressure  of  1.83  atm.  With  a 
mass  Tlow  rate  of  1.57  kg/sec,  the  inlet  velocity  and  static  temperature  were 


159  m/sec  and  543°K,  respectively;  the  inlet  Mach  number  was  0.351.  Initial 
turbulent  kinetic  energy  and  dissipation  rate  values  were  established  in  the 
same  manner  as  for  the  cold  flow  calculations,  and  the  geometry  of  the  combustor 
and  recirculation  zone  was  the  same  as  in  the  cold  flow  calculations. 

Both  combustion  efficiency  and  wall  static  pressure  distribution 
data  are  available  for  this  configuration,  and  the  results  of  the  modular 
model  calculation  of  these  quantities  are  shown  in  comparison  with  the 
experimental  data  in  Fig.  27.  The  combustion  efficiency  was  computed  from 
the  mass-average  temperature  at  each  axial  location  in  the  calculation, 
using  the  JANNAF  temperature-rise  combustion  efficiency  definition.  For 
these  data,  TT5  (ideal)  =  2178°K.  Note  that  in  these  comparisons,  data 
points  are  shown  for  the  fuel/air  ratio  0.053  value  at  which  the  computa¬ 
tions  were  carried  out.  These  data  points  were  obtained  from  plots  of  com¬ 
bustion  efficiency  vs.  fuel/air  ratio  presented  in  Ref.  26  for  three  dif¬ 
ferent  values  of  combustor  L/D.  The  band  shown  for  each  data  point  repre¬ 
sents  the  range  of  observed  combustion  efficiencies  as  a  function  of  fuel/ 
air  ratio  and  is  a  better  indication  of  the  overall  trend  of  the  combustion 
efficiency  vs.  length  data  than  are  the  individual  data  points  themselves. 

Some  caution  is  advised  in  interpreting  the  combustion  efficiency  comparison 
shown  in  Fig.  27  since  three  different  combustor  configurations  were  involved 
in  obtaining  the  data.  Thus  the  relationship  between  recirculation  zone 
length  and  combustor  length  is  different  for  each  of  the  three  combustors 
tested. 'On  the  other  hand,  the  static  pressure  data  (for  the  fuel -air 
ratios  which  bracket  the  fuel-air  ratio  used  in  the  computation)  and 
the  predictions  are  both  for  a  combustor  L/D  of  3.0. 

Overall,  the  level  of  agreement  with  the  experimental  data  is 
reasonably  good,  although  there  is  a  drop  in  both  combustion  efficiency  and 
wall  static  pressure  towards  L/D  =  3  in  the  prediction  that  is  not  evident 
in  the  data.  Both  of  these  differences  between  prediction  and  experiment 
may  be  related  to  the  inaccuracies  inherent  in  the  simple  one-step  global 
finite-rate  model  used  in  these  calculations.  Of  greater  interest,  however, 
is  the  overall  low  level  of  combustion  efficiency  shown  by  these  results, 
with  a  maximum  combustion  efficiency  of  less  than  0.6  at  L/D  =  3.  While 
this  represents  a  short  combustor,  the  premixed  inlet  condition  should 


alleviate  some  of  the  combustion  efficiency  problems  that  short,  flame- 
holderless  sudden-expansion  combustors  are  known  to  exhibit.  In  this  case 
the  computation  shows  that  the  low  combustion  efficiency  arises  from  the 
fact  that  the  combustion  is  premixed.  With  a  uniform  directed  flow  fuel/ 
air  ratio,  the  recirculation  region  is  also  at  the  same  overall  fuel/air 
equivalence  ratio  which,  in  this  case,  is  relatively  lean  ( 4»  =  0.8).  Thus 
the  recirculation  region  temperature  is  lower  than  it  would  be  at  an  equiv¬ 
alence  ratio  nearer  1.0.  Since  it  is  the  recirculation  region  that  stabi¬ 
lizes  the  flame,  and  provides  the  energy  source  for  flame  propagation  into 
the  directed  flow,  this  results  in  relatively  slow  flame  propagation  and  an 
overall  low  combustion  efficiency. 

While  these  results  indicate  a  good  overall  level  of  agreement 
with  experimental  data  for  a  premixed  combustor,  the  nonpremixed  case  in¬ 
volving  liquid  fuel  injection  is  of  considerably  greater  practical  interest. 

As  a  test  of  the  applicability  of  the  modular  approach  (which,  as  described 
in  Section  3.4,  incorporates  a  semiempirical  model  for  liquid  fuel  injection) 
to  the  analysis  of  nonpremixed  ramjet  combustor  flowfield,  the  experimental 
configuration  tested  by  Schmotolocha  and  Economos  (Ref.  27)  was  selected. 

For  these  experiments  an  area  ratio  A3/A2  =  9  combustor  was  used,  with  liquid 
JP-4  injected  perpendicularly  from  the  inlet  walls  through  an  equally-spaced 
ring  of  orifices  located  6  in.  upstream  of  the  dump  plane.  The  inlet  air 
Mach  number  for  this  combustor  was  0.49  and  the  global  fuel/air  equivalence 
ratio  was  0.63.  Since  the  experiments  involved  a  heavy  wall  combustor,  and 
the  code  is  presently  limited  to  either  adiabatic  or  isothermal  wall  boundary 
conditions,  the  effects  of  heat  loss  to  the  walls  on  the  predicted  temperature- 
rise  combustion  efficiency  were  approximated  by  using  a  lower  air  temperature 
value  measured  at  the  exit  of  the  combustor  during  isothermal  testing  in  place 
of  the  actual  inlet  air  temperature. 

Model  coefficients  were  taken  to  be  those  established  from  the 
modeling  of  the  premixed  dump  combustor  described  earlier  in  this  section, 
despite  the  large  difference  in  dump  area  ratio  between  the  combustor  used 
by  Craig,  et  al .  (Ref.  26)  and  the  present  combustor.  Results  of  these  compu- 
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tations  are  shown  in  Fig.  28.  In  the  lower  part  of  Fig.  28,  the  computed  com¬ 
bustion  temperature  distribution  is  shown.  Coalescence  of  the  isotherms  along 
the  recirculation  zone  boundary  reflects  the  highly  turbulent  shear  layer 
region.  The  predicted  equivalence  ratio  in  the  recirculation  zone,  4> RZ,  is 
1.87,  as  compared  to  a  measured  value  of  1.70,  for  the  globally  fuel-lean 
overall  equivalence  ratio  of  <J>0  =  0.63.  Note  that  the  vaporization  process 
resulting  from  liquid  fuel  injection  into  the  1110°R  airstream  effectively 
reduces  the  heterogeneous  mixture  temperature  levels,  as  shown  by  the  presence 
of  a  1000°R  isotherm.  A  comparison  of  predicted  and  measured  radial  tempera¬ 
ture  profiles  is  provided  at  the  combustor  exit  plane.  As  expected,  the 
measured  temperature  near  the  wall  is  lower  than  theoretical  due  to  local 
heat  loss.  In  the  region  near  the  centerline  the  local  temperature  is  substan¬ 
tially  underpredicted,  which  could  result  from  deficiencies  inherent  in  the 
turbulence  modeling  or  in  the  global  finite-rate  chemical  kinetics  model  (which 
substantially  overpredicts  the  required  reaction  time  compared  to  a  more  com¬ 
plete  chemical  kinetics  formulation),  or  both.  However,  the  mass  flow  in  this 
region  is  small  relative  to  the  overall  combustor  mass  flow,  so  that 
overall  combustion  efficiency,  which  is  based  on  mass-averaged  tempera¬ 
tures,  is  well  represented.  This  is  shown  in  the  upper  part  of  Fig.  28; 
the  two  experimental  points  represent  combustion  efficiency  obtained 
from  exit  species  distributions  and  from  exit  total  temperature. 

These  results  indicate  another  feature  of  the  modular  model, 
which  is  its  depiction  of  the  details  of  the  flowfield  that  result  from 
the  interaction  between  the  fuel  injection  process  and  the  distribution 
of  fuel  at  the  inlet  and  within  the  recirculation  region.  For  example, 
a  variety  of  observations  have  shown  that  with  wall  orifice  fuel  injection, 
the  recirculation  zone  equivalence  ratio  is  generally  more  fuel-rich  than 
the  overall  (global)  equivalence  ratio.  This  result  has  been  demonstrated 
by  the  data  of  Schmotolocha  and  Economos  (Ref.  27),  which  was  obtained  for 
three  different  inlet  conditions  and  at  two  different  locations  of  the  fuel 
injection  orifices  upstream  of  the  dump  plane.  Indeed,  the  data  described 
in  Ref.  27  indicate  that  for  oyerall  equivalence  ratios  greater  than  0.2, 
the  equivalence  ratio  in  the  recirculation  zone  is  always  greater  than 
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Figure  28.  Comparison  of  Modular  Model  Prediction 
with  Measured  Dump  Combustor  Perfor¬ 
mance,  Liquid  Fuel  Injection.  Data 


unity.  These  data  are  shown  in  Fig.  29,  in  which  the  results  of  modular 
model  calculations  at  two  overall  equivalence  ratios  are  also  shown. 

The  comparison  between  the  data  of  ref.  27  and  the  calculations  indicates 
good  agreement;  what  the  computational  detail  shows  is  that  this  result 
arises  from  the  fact  that  although  the  peak  equivalence  ratio  in  the  fuel- 
air  distribution  at  the  inlet  plane  is  away  from  the  wall  and  within  the 
directed  flow,  the  equivalence  ratio  near  the  wall  and  thus  adjacent  to  the 
recirculation  zone  downstream  of  the  dump  plane  is  still  substantially 
greater  than  one. 

Poppet  fuel  injectors  provide,  by  their  design,  a  fuel  penetra¬ 
tion  that  is  basically  independent  of  fuel  flow  rate  over  their  design 
range.  This  can  be  accounted  for  within  the  modular  model  approach  by 
assigning  a  fixed  value  for  fuel  penetration  distance  and  bypassing  the 
fuel  penetration  calculation  described  in  Section  2.3.4.  Fuel  mixing  and 
vaporization  are  still  accounted  for  using  the  semiempirical  approach  out¬ 
lined  in  the  latter  section.  Recent  work  carried  out  at  SAI  has  involved 
the  investigation  of  ramjet  combustor  performance  characteristics  for  en¬ 
gines  using  poppet  fuel  injectors,  for  flight  conditions  representing 
relatively  low  Mach  numbers.  Results  of  some  of  this  work  are  shown  in 
Fig.  30;  these  calculations  make  use  of  the  version  of  the  modular  model 
incorporating  the  advanced  quasiglobal  formulation.  The  comparison  shown 
in  Fig.  30  indicates  that  the  computational  results  agree  very  well  with 
the  available  data.  Since  a  fixed  radial  fuel  injection  location  was 
assumed  for  these  calculations,  the  equivalence  ratio  in  the  recirculation 
region  is  only  slightly  above  stoichiometric  and  the  temperature  is  essen¬ 
tially  the  same  in  all  cases.  The  performance  dropoff  observed  in  both 
the  experimental  data  and  the  calculational  results  as  $-*-1.0  is  a  result  of 
the  increasing  region  of  fuel-rich  equivalence  ratio  as  the  overall  (global) 
equivalence  ratio  approaches  one.  Combustion  in  this  fuel-rich  region  re¬ 
sults  in  lower  temperatures  than  are  produced  by  stoichiometric  or  slightly 
fuel-lean  combustion.  Since  the  ideal  equilibrium  temperature  is  that 
achieved  by  a  one  dimensional  combustion  process  starting  from  the  average 
combustor  inlet  temperatu)  '  and  at  the  average  fuel-air  ratio  of  the  flow, 
the  actual  mass-averaged  temperature  for  the  nonuniform  stream  will  remain 


Figure  29.  Reciriculation  Zone  Fuel  Concentration 


lower,  even  when  all  fuel  is  consumed.  Thus  there  is  a  dropoff  in  combus 
tion  efficiency  that  represents  the  nonuniformity  of  the  flow  upstream  of 
fuel  ignition  in  the  directed  flow. 

2.3.7  Discussion 

The  results  documented  in  this  section  provide  considerable 
evidence  of  the  overall  accuracy  of  the  modular  formulation  for  a  ramjet 
combustor  model.  They  also  provide  certain  insights,  particularly  with 
respect  to  the  interaction  of  the  mixing  and  combustion  processes,  that 
are  useful  in  interpreting  the  results  obtained  in  ramjet  test  programs. 
For  example,  the  results  just  discussed  indicate  that  the  high  fuel-air 
ratio  efficiency  falloff  that  is  often  seen  in  ramjet  performance  data 
results  from  the  existence  of  locally  fuel-rich  gases  when  the  overall 
fuel -air  ratio  is  stoichiometric  and  the  fuel  is  injected  just  upstream 
of  tne  combustor  inlet.  However,  the  results  for  a  premixed  combos' 
also  shown  in  this  section  show  that  at  lower  overall  fuel-air  re  a 
uniformly  mixed  inlet  stream  does  not  provide  high  levels  of  per  .,-mance 
either.  In  this  case,  because  the  inlet  flowfield  is  premixed,  am  „  an 
overall  low  equivalence  ratio,  the  flame-stabilizing  recirculation  zones 
are  also  operating  fuel-lean.  This  results  in  a  recirculation  zone  temp¬ 
erature  that  can  be  considerably  less  than  the  stoichiometric  flame  temp¬ 
erature  and  thus  a  reduced  rate  of  propagation  of  the  flame  out  of  the 
recirculation  zones  and  across  the  flow.  When  the  fuel  is  injected  from 
discrete  injectors  just  upstream  of  the  combustor  inlet,  the  modular 
model  results  show,  in  agreement  with  experiment,  that  the  recirculation 
zone  equivalence  ratio  is  essentially  always  greater  than  unity,  i.e., 
fuel-rich.  For  an  overall  fuel-lean  combustor  operation,  this  results 
in  higher  recirculation  zone  temperatures  than  would  be  achieved  in  a 
premixed  combustor  and  hence  enhanced  flame  propagation  and  improved 
combustion  efficiency. 


The  value  of  this  approach  lies  both  in  its  utility  as  a  tool 
for  the  investigation  of  the  fundamental  mechanisms  that  inderlie  observed 
ramjet  performance  trends  and  in  serving  as  a  foundation  for  further  devel¬ 
opment  of  more  sophisticated  approaches,  incorporating,  for  example,  treat¬ 
ments  of  the  two-phase  flow  in  a  spray  combustion  situation.  Indeed,  one 
of  the  considerations  involved  in  the  development  of  the  spray  flame  code 
discussed  at  the  beginning  of  this  section  is  that  it  be  compatible  with 
and  usable  as  the  directed  flow  portion  of  a  spray-combustor  modular  model. 
Combination  of  this  spray  code  and  advanced  versions  of  the  boron  particle 
combustion  treatments  with  the  modular  model  framework  outlined  in  this 
section  would  provide  a  useful  tool  for  the  investigation  of  slurry-fuel 
combustion.  Finally,  it  is  possible  to  incorporate  the  effects  of  swirl 
into  the  directed  flow  module  while  using  the  sti rred-reactor  formulation 
to  represent  the  large-scale  recirculation  regions  that  are  formed  as  a 
result  of  the  application  of  large  amounts  of  swirl  to  a  flow.  Such  an 
extension  of  the  approach  would  provide  a  useful  tool  for  the  investiga¬ 
tion  of  fundamental  phenomena  in  chemically-reacting  swirling  flowfields. 


2.3.8 


Turbulence  Models  for  Ramjet  Flowfields 


Another  aspect  of  the  ramjet  modeling  work  carried  out  under  the 
present  program  is  the  investigation  of  advanced  turbulence  models.  Although 
the  two-equation  k-e  approach  has  provided  quite  widely  useful,  and  is 
incorporated  in  the  modular  model  developed  as  part  of  this  program,  it 
does  rest  on  a  gradient  diffusion  assumption  for  momentum  transfer,  and 
questions  have  been  raised  as  to  the  appl icabi 1 i ty  of  this  approach.  The 
gradient  diffusion  assumption  is  avoided  if  the  Reynolds  stress  equations 
are  solved  directly,  but  this  adds  three  coupled  partial  differential  equa¬ 
tions  to  the  set  to  be  solved.  If  it  is  assumed,  however,  that  convection 
and  diffusion  of  each  of  the  Reynolds  stresses  can  be  related  to  convection 
and  diffusion  of  turbulence  kinetic  energy,  then  the  partial  differential 
equations  reduce  to  algebraic  expressions  and  only  the  turbulent  kinetic 
energy  (and  turbulent  kinetic  energy  dissipation  rate)  equations  remain  to 
be  solved;  their  solution  is  r.c  lo"  pi  i  shed  in  the  same  manner  as  in  the  two- 
equation  model . 
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The  assumption  that  convection  and  diffusion  of  Reynolds  stress 
can  be  related  directly  to  convection  and  diffusion  of  turbulent  kinetic 
energy  forms  the  foundation  of  the  algebraic  stress  model  (ASM)  first 
proposed  by  Rodi  (Ref.  28).  Since  each  stress  component  is  obtained  directly, 
this  approach  avoids  use  of  a  gradient  diffusion  hypothesis.  Moreover,  for 
flows  in  which  multiple  stress  components  are  apparent,  e.g.,  swirl  flows, 
the  ASM  approach  avoids  an  isotropic  eddy  viscosity  assumption  and  avoids 
any  empirical  relationship  between  different  stress  components.  This  can  re¬ 
sult  in  substantial  differences  in  overall  flowfield  prediction,  as  is 
evidenced  by  the  results  shown  in  Fig.  31.  These  results  are  for  a  constant- 
angle  swirl  flow  at  a  swirl  number  of  2,  with  identical  initial  and  boundary 
conditions  for  each  computation.  Note  that  while  both  approaches  yield  closely 
similar  outer  reciriculation  zone  lengths,  there  are  large  differences  in  the 
predicted  extent  and  location  of  the  centerline  recirculation  zones,  as  well 
as  in  the  uniformity  of  the  exit  flow.  This  latter  characteristic  is  indicated 
by  the  parallelism  of  the  stream  function  contours  at  the  exit  in  the  k-e  model 
prediction,  compared  to  the  convergence  toward  the  centerline  of  the  stream 
function  contours  predicted  by  the  ASM  model,  as  well  as  the  near-wall  value 
of  normalized  stream  function. 

Although  the  results  shown  in  Fig.  31  indicate  strong  differences 
between  the  predictions  of  the  two  turbulence  model  formulations,  there  is  as 
yet  insufficient  data  with  which  to  evaluate  the  two  approaches,  particularly 
with  respect  to  swirling  dump  combustor  flows.  However,  work  already  carried 
out  (Ref.  3)  has  shown  that  the  use  of  the  ASM  formulation  does  improve  pre¬ 
dictions  in  planar  recirculating  flows  with  large  amounts  of  streamline  curva¬ 
ture,  and  indirect  evidence  (such  as  observations  of  recirculation  zone  size 
and  levels  of  exit-plane  nonuniformity)  also  indicates  that  further  develop¬ 
ment  of  the  ASM  formulation  should  be  pursued. 

A  second  aspect  of  the  comparison  of  turbulence  model  results  for 
swirling  ramjet  combustor  flowfields  involves  the  examination  of  the  effects 
of  different  initial  swirl  distributions.  An  example  of  this  is  provided 
by  the  results  shown  on  Figs.  32  and  33,  which  show  radial  profiles  of  axial 
(Fig.  32)  and  tangential  (Fig.  33)  velocity,  for  a  swirling  flow  in  a  1.5:1 
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33.  Tangeriial  Velocity  Profiles,  S  =  0.4, 
Different  Initial  Swirl  Distributions 


diameter  ratio  sudden  expansion  for  four  different  types  of  inlet  swirl. 

Note  that  the  low  area  ratio  was  chosen  to  minimize  inertial  effects  and 
thus  maximize  the  effects  of  the  turbulence  model  formulation  on  the  flowfield 
prediction.  The  overall  swirl  number  in  this  demonstration  is  held  constant 
at  0.4,  and  there  is  a  centerbody  at  the  origin  of  this  flow  (within  the 
inlet).  A  nonuniform  37  (axial)  by  26  (radial)  grid  was  used,  and  predictions 
are  shown  for  both  a  two-equation  and  ASM  turbulence  model.  Despite  the 
small  size  of  the  figures,  it  is  immediately  apparent  from  the  tangential 
velocity  profiles  shown  in  Fig.  33  that  the  ASM  approach  prediction  of  the 
downstream  mixing-out  process  differs  substantially  from  that  predicted  using 
the  two-equation  approach.  It  was  noted  earlier  that  there  is  a  tendency  in 
the  work  reported  in  the  literature  for  predictions  of  swirl  flows  to  relax 
to  a  profile  characteristic  of  a  forced  vortex,  w«r,  whatever  the  initial 
tangential  velocity  distribution.  Such  a  tendency  is  observable  in  the 
profiles  predicted  using  the  k-e  approach  shown  on  Fig.  33(a),  but  the  ASM 
predictions.  Fig.  33(b)  do  not  exhibit  thi^  tendency.  This  effect  thus 
appears  to  be  directly  related  to  the  existence  of  an  isotropic  eddy 
viscosity  assumption  implicitly  in  the  k-e  model  application,  whereas  with 
the  ASM  (different)  effective  viscosities  are  obtained  for  each  of  the  three 
shear  stress  components. 

2.4  MODELING  OF  THE  DUCTED  ROCKET  COMBUSTION  PROCESS 

The  work  being  carried  out  in  this  program  under  the  task  area 
outlined  in  Section  1(d)  had  as  its  objective  the  development  of  a  new 
theoretical  model  to  describe  the  three-dimensional  reacting  flow  typical 
of  multiple-inlet  ducted-rocket  combustion  chambers.  It  was  supported, 
through  AFOSR,  by  the  Ramjet  Division  (RJT),  AFWAL.  Because  of  the  paucity 
of  experimental  data  on  the  aerodynamics  of  the  ducted  rocket  configuration, 
the  modular  model  development  work  under  this  task  area  utilized  the 
hierarchy  of  models  concept,  in  which  unified  elliptic  aerodynamic  models 
were  used  to  define  characteristic  flowfield  regions  for  the  development  of 
modular  model s  which  incorporate  the  fundamental  chemical  kinetic  formula¬ 
tions  required  to  address  the  performance  aspects  of  ducted  rockets. 


Aerodynamic  calculations  were  carried  out  for  a  three-dimensional 
configuration  similar  to  an  experimental  ducted  rocket  combustor  tested 
at  AFWAL/RJT.  The  three-dimensional  elliptic  analysis  used  for  these  com¬ 
putations  was  a  modified  version  of  an  analysis  similar  to  the  Imperial 
College  TEACH  code,  incorporating  a  two-equation  turbulent  kinetic  energy 
turbulence  model  similar  to  that  used  for  the  computations  described  in 
Sections  2.1  and  2.3.  A  three-dimensional  grid  is,  of  course,  required, 
with  21  axial,  10  radial,  and  18  circumferential  grid  points.  This  grid  is 
not  fine  enough  to  provide  grid-independent  results,  but  since  it  involves 
3780  node  points  it  is  a  compromise  between  the  fineness  needed  to  obtain 
grid-independence  and  the  limitations  of  computer  storage  and  speed*.  The 
configuration  examined  involves  two  symmetrically-disposed  inlets  at  180° 
apart:  this  is  not  identical  to  the  AFWAL  configuration  but  was  selected 
to  simplify  the  implementation  of  boundary  conditions  in  these  computations. 

A  sketch  of  this  combustor  geometry  is  shown  in  Figure  34,  from  which  it  can 
further  be  seen  that  the  inlets  are  disposed  at  an  angle  of  45°  to  the 
longitudinal  axis  of  the  combustor.  The  air  inlet  velocity  was  selected 
to  correspond  to  an  airflow  rate  of  2  lb/sec/inlet. 

The  velocity  vectors  computed  for  this  flow  field  in  a  plane 
through  the  inlets  and  along  the  combustion  centerline  are  shown  in  Fig.  35. 
Two  recirculation  zones  can  be  seen,  one  at  the  head  end  of  the  combustor 
and  the  second  along  the  outer  walls  just  downstream  of  the  inlets.  Contours 
outlining  the  approximate  sizes  of  these  recirculation  regions  are  shown  on 
Fig.  35.  For  fuel  injection  in  the  inlet  arms  (i.e.,  the  side-dump  liquid 
fueled  ramjet  configuration)  both  of  these  recirculation  regions  are  poten¬ 
tial  flameholding  sites.  However,  if  fuel  is  injected  from  the  combustor 
head  end  (the  ducted  rocket  configuration)  it  would  be  expected  that  the 
primary  flameholding  region  would  involve  the  recirculation  zone  at  the 
combustor  head  end. 


*The  solution  procedure  involves  eight  variables:  three  velocity  components, 
the  pressure,  the  flow  enthalpy,  and  a  species  used  as  a  tracer,  as  well  as 
the  turbulent  kinetic  energy  and  its  dissipation  rate.  Thus  for  the  grid 
noted,  30,240  locations  are  required  for  the  storage  of  the  dependent  var¬ 
iables  alone,  for  this  nonreacting  (but  possibly  nonisothermal )  calculation. 


The  radial  extent  of  the  recirculation  zone  at  the  head  end  of  the 
ducted  rocket  combustor  is  depicted  in  Fig.  36.  Of  interest  in  this  view 
are  the  distinctive  four-lobed  shape  of  the  recirculation  region  and  the 
stagnation  point  shown  on  the  flowfield  centerline.  Through  computations 
such  as  represented  by  Figs.  35  and  36  it  is  possible  to  build  up  a  three- 
dimensional  view  of  the  size  and  shape  of  the  recirculation  zones  in  the 
ducted  rocket  combustor.  This  provides  one  of  the  required  pieces  of  infor¬ 
mation  for  the  development  of  a  modular  model  of  this  flowflield. 

Various  types  of  modular  models  can  be  postulated  for  a  flowfield 
such  as  shown  in  Figs.  35  and  36.  However,  one  of  the  most  useful,  because 
of  its  computational  speed  and  flexibility,  is  the  well-stirred  reactor/plug 
flow  with  distributed  secondary  injection/entrainment  approach.  In  this  model 
the  intense  mixing  that  occurs  in  the  region  of  the  inlets  in  the  ducted 
rocket  combustor  is  modeled  using  a  well-stirred  reactor  formalism,  while 
the  subsequent  mixing  and  combustion  of  fuel  and  air  not  transported  into  the 
recirculation  zone  is  modeled  as  a  plug  flow  with  distributed  entrainment. 
While  this  approach  is  fast  and  flexible,  it  requires,  in  addition  to  a 
specification  of  the  size  and  shape  of  the  recirculation  zone,  a  measure 
of  the  proportion  of  the  flow  involved  in  the  recirculation  region  and  in 
the  subsequent  mixing  region.  This  specification  of  the  proportion  of  the 
flow  found  in  different  regions  can  be  provided  using  the  elliptic  aerody¬ 
namic  formulation  with  a  species  tracer. 

Calculations  carried  out  to  date  to  examine  the  mixing  and  kinetics 
limitations  within  the  side-dump  combustor  used  at  AFWAL  have  not  been 
successful.  In  part  this  state  of  affairs  arises  from  the  complexity  of  the 
overall  flowfield,  but  in  general  the  characteristics  of  the  combustor  appear 
to  be  such  that  kinetics  are  not  limiting.  Mixing-limited  situations  are 
inherently  more  difficult  to  approach  with  the  sort  of  modular  model  being 
applied  in  this  phase  of  the  program.  The  work  in  this  specific  area  is 
continuing  under  a  follow-on  AFOSR  contract. 

It  should  be  noted  that  the  configuration  tested  at  AFWAL  is  not 
a  true  ducted  rocket  but  represents  a  side-dump  ramjet  combustor  configura¬ 
tion.  One  of  the  essential  features  of  a  ducted  rocket  is  the  presence  at 
the  head  end  of  the  combustor  of  one  or  more  jets  of  higher  temperature. 


Figure  35.  Longitudinal  Cross-Section  through 
Inlet  Arms,  Ducted  Rocket  Combustor 
Aerodynamic  Flowfield 


View  B-B  (See  f igu re  34) 


Figure  36.  Lateral  Cross-Section  through  Inlet 
Arms,  Ducted  Rocket  Combustor  Aerody 
namic  Flowfield 


fuel-rich,  partially  reacted  products  of  the  rocket  combustion  process. 

These  jets  are  sonic,  or  more  usually,  highly  underexpanded  supersonic  flows. 
Thus  the  flowfield  within  a  true  ducted  rocket  involves  both  subsonic  and 
supersonic  regions,  with  strong  shock  and  pressure  wave  interactions,  and 
could  be  expected  to  be  considerably  different  in  character  from  the  flow- 
field  shown  in  Figs.  35  and  36.  This  aspect  has  not  been  accounted  for  in 
any  ducted  rocket  modeling  so  far  described  and  will  be  considered  in  con¬ 
tinuing  SAI  work  in  this  area. 
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out  under  the  task  area  outlined  in  Section  1(d)  involves  close  interaction 
with  experimental  work  carried  out  by  Drs.‘F.D.  Stull  and  R.D.  Craig  at  AFWAL/ 
ROT.  In  addition,  the  boron  slurry  modeling  effort  described  in  Section  1(c) 
involves  coordination  with  experimental  work  being  done  by  Dr.  Klaus  Schadow 
at  NWC. 

In  addition.  Dr.  Edelman  served  on  the  Propuslion  Technology  Panel  of 
the  Military  Space  System  Technology  Model  (MSSTM)  evaluation  activity  for 
the  Space  Division  of  the  Air  Force  Systems  Command,  and  Drs.  Edelman  and 
Harsha  have  both  been  active  in  participating  in  JANNAF-sponsored  workshops 
on  swirl  flows  in  ramjets,  two-phase  flow  in  ramjet  combustors,  and  ducted 
rocket  combustion  technology. 


FILMED 


1-85 

9 


DTIC 


